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ABSTRACT 


A  micromechanical  model  is  advanced  in  order  to  study  the  stress  transfer  and 
associated  damage  and  failure  in  classes  of  conventional  and  textile  type  fibrous 
composites.  Unidirectionally  reinforced  matrix  with  straight  and  undulated  fibers  define 
the  repeating  constructing  cell  for  conventional  and  textile  composites,  respectively. 
Starting  with  the  case  of  straight  reinforcement,  we  approximate  and  model  the  actual 
discrete  composite  as  a  concentric  cylindrical  system.  For  axisymmetric  loading,  and 
upon  adopting  some  appropriate  restrictions  on  the  radial  behavior  of  some  field 
quantities,  an  elasticity-based  procedure  reduces  the  two-dimensional  field  equations, 
which  hold  in  both  fiber  and  matrix  components  together  with  the  appropriate  interface, 
symmetry  and  boundary  conditions,  to  a  quasi-one-dimensional  system.  This  analysis  is 
further  extended  to  cases  involving  undulated  fibers.  Based  upon  local  directions  (slopes) 
of  the  undulated  fibers,  the  linear  transformation  is  used  to  obtain  local  stress 
distributions  along  the  undulated  fibers.  The  total  stress  field  is  found  to  be  combinations 
of  these  local  stresses  and  the  inherent  contributions  obtained  from  the  transformations  of 
the  normal  loads  along  the  undulated  directions  in  the  absence  of  reinforcement.  This 
simple  system  retains  total  account  of  the  system’s  physics  and  presents  itself  in  the  form 
of  coupled  partial  differential  equations  in  the  longitudinal  displacements  and  stresses  of 
both  the  fiber  and  matrix  components. 

According  to  this  model,  damage  is  simulated  in  the  form  of  stress  free  boundary 
conditions.  Perpetuation  of  damage  is  based  upon  the  maximum  normal  stress  criterion. 
The  adverse  effect  of  such  damage  on  the  stiffness  properties  of  the  composite  is 


predicted.  Results  show  the  favorable  effect  of  undulation  in  decreasing  the  rate  of 
property  degradation  with  increasing  damage.  The  model  is  quite  general  and  has  been 
applied  to  several  situations.  These  include  response  to  static  loading,  and  other  vibration 
and  wave  propagation  applications.  Confidence  in  our  modeling  procedure  has  been 
confirmed  by  comparisons  with  whatever  available  analytical  models,  experimental  data 


and  finite  element  calculations. 
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CHAPTER  1 


INTRODUCTION 

Textile  composites  are  emerging  as  a  new  generation  of  advanced  complex 
materials  with  the  potential  of  enhancing  shortcomings  of  conventional  unidirectional 
fibrous  composites.  Advanced  structural  applications  have  recently  required  significant 
improvements  of  the  damage  tolerance  and  reliability  of  composites.  To  achieve  such 
improvements,  a  high  level  of  through  thickness  and  interlaminar  strengths  is  required. 
Such  requirements  can  be  easily  tailored  by  undulated  (textile)  fibrous  composites.  For 
these  classes  of  complex  materials,  only  limited  and  mostly  recent  attempts  have  been 
carried  out  on  their  modeling  and  analysis.  In  micromechanical  modeling,  as  a  minimum, 
models  are  needed  to  characterize  their  effective  properties.  Even  here,  very  little  work  is 
available.  For  more  sophisticated  and  advanced  applications,  especially  those  concerned 
with  damage  tolerance,  information  about  the  local  distributions  of  stresses  in  the  fiber 
and  matrix  constituents  individually  are  practically  absent  although  of  prime  importance. 
To  put  future  discussion  into  perspective,  it  is  instructive  at  this  stage  to  digress  and  give 
a  short  account  of  the  various  classes  of  composite  materials.  Composite  materials  have 
gone  several  classification  processes,  both  on  the  micro-  and  on  the  macro-scales.  One 
classification  in  particular  is  based  on  the  shape  and  distribution  of  one  phase,  namely  a 
fiber,  into  another  phase,  the  matrix.  In  this  context,  fibrous  composites  are  divided  into 
straight  fibrous  composites  and  undulated  fibrous  composites.  This  latter  class  can  exist 
in  several  shapes,  including  2-D  and  3-D  braids,  plain  weave  and  3-D  woven  composites, 
all  define  classes  of  textile  materials. 
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Textile  structural  composites  in  particular  are  defined  as  textile-containing 
materials  designed  for  structural  or  load  bearing  applications.  Like  conventional 
composites,  because  of  their  higher  strength-to-weight  ratio  and  stifFness-to-weight  ratio, 
textiles  are  also  considered  superior  to  ordinary  metal  alloys.  The  idea  of  using  an 
assembly  of  yams  woven  together  to  produce  a  stmctural  composite  came  from  the 
clothing  industry.  Even  the  technologies  used  in  producing  such  composites  have  their 
origins  in  clothing  factories.  As  textiles  continue  to  be  more  and  more  important  in  many 
stmctural  applications,  there  is  a  growing  need  for  accurate  prediction  of  their  mechanical 
performance  under  various  loading  conditions. 

It  is  usually  the  case  that  analysis  of  textile  composites  is  based  on  analyses  of 
straight  fibrous  composites.  Comparatively  speaking,  for  composites  reinforced  by 
straight  fibers,  the  literature  is  almost  saturated  by  the  number  of  studies  that  are  devoted 
to  deriving  their  effective  properties.  It  is  beyond  this  work,  however,  to  review  this  vast 
literature.  For  updated  review  see,  for  example,  Nayfeh  [1]  and  Al-Huniti  [2].  Studies 
concerning  the  distribution  of  local  stresses  and  the  degree  of  load  transfer  between  the 
straight  fiber  and  matrix  components  have  also  received  reasonable  attention.  The 
importance  of  the  problem  is  attributed  to  the  fact  that  accurate  prediction  of  the  resulting 
stress  distribution  in  both  constituents  has  a  strong  bearing  upon  the  predicted  mechanical 
integrity  of  the  loaded  composite. 

The  majority  of  existing  models  of  straight  fibrous  composites  included  various 
restrictions  on  the  material  properties,  field  equations  and  interface  and  boundary 
conditions.  According  to  some  theoretical  models,  these  resulted  in  a  second  order 
differential  equation  that  describes  the  behavior  of  either  interfacial  shear  stress  or  normal 
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stress  in  the  fiber  or  the  matrix;  examples  include  the  work  of  Smith  and  Spencer  [3], 
Nayfeh  [4,5],  Steif  [6],  Gao  et  al.  [7],  Carapella  et  al.  [8],  and  Kishore  et  al.  [9], 

More  recently,  McCartney  [10]  presented  a  micromechanically-based  analysis 
that  can  predict  the  stress  transfer  between  the  fiber  and  matrix  in  a  uniaxially  straight 
fiber  reinforced  composite.  The  system  was  modeled  as  a  concentric  cylindrical  system 
with  stress  free  outer  boundary.  Broadly  speaking,  this  situation  corresponds  to  the  case 
of  plane  stress  loading.  Thermomechanical  loads  in  which  the  mechanical  ones  are 
applied  as  tension  along  the  fiber  direction  are  applied.  The  approach  uses  an  elasticity 
based  analysis  and  starts  with  approximating  some  of  the  field  variables  and  proceeds  to 
satisfy  the  appropriate  field  equations  of  equilibrium,  constitutive  relations  and  the 
associated  interface  and  boundary  conditions.  In  contrast  to  above  mentioned  older 
models,  McCartney’s  analysis  lead  to  a  fourth  order  differential  equation  that  has  more 
versatility  leading  to  better  prediction  of  the  actual  stress  distribution  in  the  composite.  In 
particular,  McCartney  was  able  to  simulate  localized  damage  by  stress  free  boundary 
conditions. 

An  even  fewer  number  of  attempts  are  available  which  model  the 
micromechanical  distribution  of  stresses  in  other  than  straight  fiber  reinforced 
composites;  namely  in  classes  textile  materials.  Due  to  the  complexity  of  the  system, 
some  finite  element  and  boundary  element  based  studies  have  been  applied  to  the 
problem  at  hand  to  obtain  approximate  numerical  solutions.  Nedele  and  Winsom  [11],  in 
a  finite  element  parametric  study,  showed  that  geometrical  variations  in  unidirectional 
composites  have  a  higher  influence  on  the  stresses  than  the  changes  in  the  constituent 
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properties.  Goldberg  and  Hopkins  [12]  developed  a  new  boundary  element  formulation 
for  the  micromechanical  analysis  of  composites  and  applied  it  to  some  textile  composites. 

Based  upon  the  above  various  models,  several  authors  have  constructed  similar 
models  for  laminated  composites  and  then  proceeded  to  estimate  the  strength  of  the 
overall  composite.  This  greater  body  of  research  has  been  devoted  to  laminated 
composites  since  they  were  the  first  to  be  used  in  industry  and  because  of  their  relative 
simplicity  of  modeling.  Failure  criteria  that  are  based  on  the  maximum  stress,  the 
maximum  strain  or  the  more  general  Tsai-Wu  theory  [13]  have  been  employed.  Among 
the  earliest  models  stands  the  work  of  Reifsnider  and  Talug  [14]  in  which  they  derived  a 
simple  shear  lag  model  to  predict  failure  in  composite  laminates.  Later  on  Highsmith  and 
Reifsnider  [15]  verified  these  results  experimentally.  More  thorough  analyses  followed, 
including  the  work  of  Hashin  [16]  and  Talreja  [17,18].  Hashin  took  the  analysis  a  step 
further  by  dealing  with  the  three  dimensional  problem  of  orthogonally  cracked  laminates 
[19],  while  the  modeling  of  McCartney  [20]  enabled  him  to  study  the  stress  transfer 
process  in  repetitively  damaged  composites  on  the  micro-level.  Damage  accumulation 
studies  in  laminated  composites  was  discussed  by  Akshantala  and  Talreja  [21]. 

A  very  limited  number  of  publications  combine  the  geometric  and  analytic 
complexity  of  straight  or  undulated  fibers  in  attempts  to  model  damage  accumulation. 
Among  these  stands  the  work  of  Allen  and  Harris  [22]  who  developed  a  model  to  predict 
the  degradation  in  the  material  properties  of  damaged  textile  composites  subject  to 
thermomechanical  loading.  Later  on,  Pochiraju  et  al.  [23]  used  the  concept  of 
parameterized  Macro-cells  to  represent  the  woven  fiber  architecture  and  then  adopted  an 
averaging  technique  to  estimate  the  textile  composite  strength  under  in-plane  loading.  To 
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our  knowledge,  there  exist  no  further  studies  of  the  damage  accumulation  or  property 
degradation  in  these  classes  composites. 

1.1  Building  Block  Approach 

The  main  objective  of  this  work  is  to  develop  a  unified  approach  to  model  the 
micromechanical  stress  distribution  in  straight  and  undulated  fibrous  composites 
subjected  to  thermal  and  mechanical  loading.  This  will  be  sought  in  the  absence  or 
presence  of  damage  mechanisms.  Damage  perpetuation  and  accumulation  leading  to 
failure  is  also  among  our  prime  objectives.  To  accomplish  these  objectives,  an  elasticity 
based  procedure,  that  involves  use  of  the  interface  quantities  and  averaging  technique,  is 
used  to  reduce  the  two-dimensional  dynamic  field  equations  together  with  the  interface, 
symmetry  and  boundary  conditions  to  a  quasi-one-dimensional  system.  Following  a 
building  block  approach,  we  specifically  consider  the  following  steps  towards  our 
understanding  of  the  undulated  system’s  behavior; 

1 .  construct  and  verify  a  model  for  the  micromechanical  stress  distribution  in  straight 
fibrous  composites  with  or  without  local  damage, 

2.  construct  and  verify  a  model  for  the  micromechanical  stress  distribution  in  undulated 
fibrous  composites  with  or  without  local  damage, 

3.  use  the  above  results  to  study  the  failure  behavior  of  these  classes  of  composites 
under  static  and  dynamic  loading. 

In  chapters  2  and  3,  we  lay  the  foundation  of  the  work  by  starting  with  the  most 
commonly  used  fiber-matrix  arrangement,  namely  the  straight  fibrous  composites.  The 
composite  repeating  unit  cell  is  modeled  as  a  concentric  cylindrical  system.  To  maintain 
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material  generality,  the  fiber  and  the  matrix  will  be  allowed  to  have  anisotropic  or 
isotropic  properties.  Two  structural  situations  can  be  dealt  with  using  the  composite 
cylinder  model.  These  relate  to  different  constraining  boundary  conditions.  On  a  macro¬ 
scale,  a  composite  panel  can  be  either  thick  or  thin.  Both  situations  correspond  to  a 
certain  lateral  constraint  condition  on  the  corresponding  unit  cell.  In  both  situations,  the 
composite  might  be  damage  free  or  has  a  single  or  multiple  fiber  breaks  or  matrix 
transverse  cracks,  or  even  a  delamination  region  between  the  fiber  and  the  matrix 
constituents.  Furthermore,  any  of  the  above  mentioned  arrangements  might  be  subjected 
to  either  a  static  or  a  dynamic  (harmonic)  type  of  loading.  These  different  possibilities  are 
arranged  below  in  figure  1.1. 

Chapter  2  is  devoted  to  the  general  composite  described  by  the  bold  line  in  the 
figure,  namely,  the  laterally  constrained  damage-free  composite  under  dynamic  loading. 
In  chapter  3,  we  exploit  the  generality  of  the  above  model  and  show  how  solutions  for 
many  other  cases  can  be  obtained  from  this  analysis.  We  then  present  comparisons  of 
these  solutions  with  published  results. 

In  chapter  4,  the  additional  geometric  complexity  of  undulation  is  added  to  the 
model  obtained  in  chapter  2.  This  will  be  done  as  follows:  for  a  given  geometric 
undulation,  we  identify  and  analjdically  describe  the  local  tangent  to  the  fiber. 
Subsequently,  we  transform  the  global  coordinates  and  applied  loads  to  the  local  systems 
consisting  of  the  tangent  directions  and  their  normal  in  the  planes  of  undulation.  We  then 
adapt  the  results  obtained  for  straight  fibers  to  “straight”  segmented  fibers  along  the 
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Straight  flbrous  composites 


Isotropic  fiber  Anisotropic 

and  matrix  fiber  or  matrix 


Laterally  free  Laterally  constrained 

composite  composite 


A  composite  with  A  composite  with  A  composite  with  Damage-free 
broken  fibers  cracked  matrix  imperfect  bonding  composite 

_ I _ 

r  I 

Static  loading  Dynamic  loading 


Figure  1.1:  Classification  of  possible  straight  fibrous  composite  arrangements. 
Bold  case  indicates  model  of  chapter  2 


tangents.  The  local  stresses  obtained  using  this  procedure  are  then  supplemented  with  the 
stress  contributions  that  inherently  rise  in  any  oriented  direction  with  respect  to  the 
loading  direction.  As  in  chapter  2,  the  analysis  of  the  undulated  fibrous  case  will  be 
carried  out  for  a  damage  free  composite  under  dynamic  loading.  In  chapter  5,  this 
analysis  will  then  be  utilized  to  solve  different  problems  involving  different  types  of 
damage,  see  fig  1.1,  under  both  types  of  loading.  Comparisons  of  the  results  with  finite 
element  solutions  will  also  be  established. 

In  chapter  6,  we  carry  the  modeling  further  and  use  the  stress  tensor  components, 
as  predicted  in  chapters  2  and  4  in  order  to  investigate  the  failure  behavior  of  straight  and 
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undulated  fibrous  reinforced  panel.  In  particular,  the  degradation  of  the  composite’s  in¬ 
plane  properties,  represented  by  Young’s  modulus  of  elasticity,  is  investigated.  Load- 
damage  dependence  is  investigated  to  determine  the  strength  reduction  of  the  composite. 
Attention  is  also  given  to  the  problem  of  delamination  progression  between  the  fiber  and 
the  matrix  and  its  effect  on  the  overall  composite  stiffness. 

Finally,  in  Chapter  7,  we  show  the  versatility  of  the  suggested  formulation  by 
applying  it  to  investigate  the  vibration  characteristics  of  damaged  fibrous  composites. 
The  effects  of  damage  location  and  the  fiber  volume  fraction  on  the  composite  system 
natural  frequencies  are  studied,  and  several  interesting  observations  are  drawn.  We  then 
conclude  by  presenting  conclusions  and  recommendations  for  further  studies. 
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CHAPTER  2 


DEVELOPMENT  OF  THE  MODEL  FOR  STRAIGHT  FIBROUS 
COMPOSITES 

2.1  Introduction 

It  was  established  in  Chapter  1  that  the  basic  building  block  for  the  modeling  of 
textile  fibrous  composites  starts  with  straight  fiber  reinforcement.  Now,  we  start  by 
developing  a  unified  generic  but  improved  analytical  model  for  the  mechanically  induced 
stress  and  displacement  fields  in  a  unidirectional  straight  fibrous  composite  under  inplane 
dynamic  loading.  In  chapter  3,  we  then  discuss  applications  of  this  model  and  verify  its 
validity  by  comparisons  to  published  analytical  results  and  finite  element  calculations. 

Several  existing  models  have  been  used  to  describe  the  behavior  of  a  fiber-matrix 
composite  mixture  under  inplane  static  loading,  see  for  example  references  [5]  and  [6]. 
The  major  improvement  in  the  present  modeling  over  existing  ones  is  the  retaining  of  the 
term  representing  the  rate  of  change  of  the  radial  displacement  with  the  axial  coordinate 
in  the  shear  constitutive  relation.  Furthermore,  since  in  practical  situations,  the  fiber  may 
be  multi-layered,  the  continuum  mixture  treatment  has  to  be  applicable  to  anisotropic 
media.  For  this  reason,  we  shall  treat  both  the  fiber  and  matrix  constituents  as 
transversely  isotropic.  Situations  involving  isotropic  material  fibers  or  matrix  can  be 
easily  obtained  as  special  cases  by  merely  imposing  the  appropriate  restrictions  on  their 
anisotropic  properties.  The  analysis  will  proceed  first  assuming  perfectly  bonded 
constituents  with  no  damage.  Later  on,  different  types  of  composite  damage  will  be 
addressed.  These  will  include  breakage  of  the  fibers  in  planes  perpendicular  to  the 
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composite  principal  direction,  transverse  cracking  in  the  matrix  and  frictional  slip  on  the 
fiber-matrix  interfaces. 

2.2  Formulation  of  the  Problem 

Consider  the  general  unidirectionally  reinforced  periodic  composite  shown  in 
figure  2.1a.  The  fiber  can  consist  of  several  concentric  layers,  and  an  additional  outer 
layer  works  as  an  interphase.  It  is  due  to  failure  in  this  layer  that  fiber-matrix  debonding 
occurs  and  frictional  slip  takes  place.  The  composite  has  a  hexagonal  symmetry  with  a 
representative  unit  cell  as  shown  in  figure  2.1b.  The  specific  geometry  shown 
corresponds  to  a  widely  common  sample.  For  our  present  analysis,  the  fiber  layers  are 
assumed  to  be  perfectly  bonded  together  and,  in  turn,  to  the  matrix.  Since  the  continuum 
mixture  solutions  are  most  appropriate  for  two  homogenous  phase  media,  namely  for 
media  consisting  of  single  material  fiber  and  matrix,  we  need  to  model  the  multilayered 
fiber  as  an  effective  homogenized  one.  Due  to  the  biased  geometry  of  the  fiber  layers,  the 
effective  fiber  will  have  anisotropic  properties  or  more  precisely,  will  exhibit  transverse 
isotropy.  The  effective  anisotropic  properties  of  the  three-component  fiber  are 
constmcted  by  a  repetitive  use  of  the  procedure  developed  by  Nayfeh  [1]  for  deriving 
effective  properties  for  two  component  fibrous  composites.  To  facilitate  our  analysis,  we 
further  approximate  the  unit  cell  with  a  concentric  cylinder  as  shown  in  figure  2.1c.  The 
unit  cell  has  an  outer  radius  with  the  effective  fiber  having  a  radius  r, . 
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Fig.  2. 1 :  The  composite  panel  model:  (a)  the  composite  panel, 

(b)  the  representative  unit  cell,  (c)  the  effective  fiber  concentric  cylinder  model. 


For  the  present  modeling,  the  behavior  of  the  composite  under  dynamic  loading  is 


described  by  the  following  two-dimensional  momentum  equations: 


da,  15,  ,  d^u 
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(2.2) 


Here  a^,  a^,  tr^and  <t„  are  the  non-vanishing  components  of  the  stress  tensor;  u  and  v 
are  the  displacement  components  in  the  longitudinal  and  radial  directions,  respectively. 
Identifying  the  fiber  and  matrix  materials  with  the  subscripts  /  and  m,  respectively,  the 
fiber  and  matrix  have  the  transversely  isotropic  materials  constitutive  relations  take  the 
form: 
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are  the  axisymmetric  strain  displacement  relations.  We  chose  to  expand  the  above 
equations  for  further  use.  For  axisymmetric  dynamic  behavior,  the  applicable  constitutive 
relations  are 
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which  hold  for  the  fiber  if  C„  =  f„  and  for  the  matrix  if  C,  =  w,. .  The  use  of  the  stress- 
strain  relations  is  more  convenient  in  the  following  analysis  than  the  commonly  used 
choice  of  strain-stress  relations. 

Several  symmetry  and  continuity  conditions  have  to  be  satisfied  in  the  model.  The 
composite  symmetry  requires  the  vanishing  of  the  shear  stresses  at  the  center  and  outer 
surface  of  the  unit  cell.  Also,  the  lateral  constraint  assumption  implies  the  vanishing  of 
the  radial  displacement  at  its  outer  surface.  Finally,  axial  and  radial  stresses  and 
displacements  have  to  be  continuous  across  the  fiber-matrix  interface.  The  above  field 
equations  are,  therefore,  supplemented  with  the  following  conditions  (see  figure  2.2): 
Symmetry: 

v^(z,0,0  =  0  ,  a-,,^(z,0,0  =  0  ,  =  0  (2.10a) 

lateral  constraint: 

'’.(z.'i,0=0  (2,10b) 
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Continuity  (on  r  =  ): 

Or^(z,r^,t)  =  o^iz,ri,t) 

o^(z,r^,t)  =  o^{z,r^,t) 


V/(z,r,,0  =  v„(z,r,,0,  (2.10c) 

Uf(z,r^,t)  =  u„(z,r^,t).  (2.10d) 


Fig  2.2:  Symmetry  and  continuity  condition  in  the  representative  unit  cell. 


2.3  Method  of  Solution 

Following  the  procedure  outlined  in  [4],  if  equations  (2.1)  and  (2.6)  are  averaged 
according  to 
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and  if  the  symmetry  and  continuity  conditions  (2.10a,  b,  c)  are  used,  one  obtains  the 
quasi-one-dimensional  momentum  equations 
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and  the  constitutive  relations 
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Here  w.  and  m_  stand  for  the  volume  fractions  of  the  fiber  and  matrix  materials, 
respectively  and  t  and  S  denote  the  momentum  and  the  constitutive  relations  interaction 
terms,  respectively;  these  quantities  are  defined  as: 

«/  =  'i"  / '■2"  >  Wm  =  -  'i") /  ^2^ ,  (2. 16a,  b) 

T  -  In  fa*  /  r,  S  =  In^v*  /  r,  (2. 17a,  b) 


where  <j*and  v*are  the  interface  shear  stress  and  radial  displacement;  namely: 

a*  =  a„f{z,  r„  /)  =  a^{z,  /■„/),  (2.18) 

V*  =  v^(z,r,,/)  =  v„(z,/i,0.  (2.19) 

The  remaining  constitutive  relations  (Eqs.  2.7-2.9)  will  now  be  utilized  to  derive 
the  relevant  partial  differential  equations  describing  the  system  in  terms  of  the  interaction 
terms  r  and  S  and  the  axial  displacements  Uf  and  i7„ .  Once  this  is  done,  the  problem 
described  by  Eqs.  (2.12-2.15)  reduces  to  a  quasi-one-dimensional  one,  which  retains  the 
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integrity  of  the  propagation  in  the  individual  fiber  and  matrix  components  subject  to  the 
interaction  terms  rand  S.  This  cannot,  however,  be  done  without  adopting  appropriate 
approximations.  These  concern  the  radial  dependence  of  the  radial  displacement  and  the 
shear  stress,  which  can  be  adopted  from  known  elasticity  solutions  (see  for  example  [4]) 


as: 

v^(z,r,t)=B(z,t)r,  (2.20a) 

=  A(z,  t)r ,  (2.20b) 

v.(z,r.()  =  B(z.O— (4-1)'',  (2.21a) 

m 

n„  r 


chosen  to  automatically  satisfy  the  symmetry  conditions  (2.10a,  b)  and  the  interface 
continuity  conditions  (2.10c). 

Introducing  these  approximations  into  the  functions  defined  in  Eq.  (2.17)  yields 


r  =  IrijA , 

(2.22) 

S  =  2n^B . 

(2.23) 

Using  these  expressions  in  equations  (2.12-2.15)  leads  to 

da,f  d^Uf 

”,  g/=-2V  +  »/P/  . 

(2.24) 

n  =  2nfA  +  n  p  ^  , 

(2.25) 

_ 

^f^zf  ““  5 

(2.26) 
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(2.27) 


aw. 


=  -2m„nfB  +  n„m,, 

m  zm  J  mil 


At  once,  we  notice  that  two  coupled  equations  can  be  obtained  by  substituting  <j^f  and 
from  Eqs.  (2.26)  and  (2.27)  into  Eqs.  (2.24)  and  (2.25),  respectively.  This  yields 

(2.28) 


-  55  .  .  a^My 


a^«m  -  55  _  .  5^m„ 

«™Wii  — ^~2w./Wi,  —  =  2nfA  +  n„p„  — ^ 

”  "  5z^  5z  ^  5/^ 


(2.29) 


Next,  we  start  by  satisfying  the  shear  constitutive  relation  (2.9).  Multiplying  it 


with  and  integrating  by  parts  in  accordance  with  Eq.  (2. 1  la),  we  directly  obtain 

.  _  ^  a5, 

u  -Ur=  —  {- — — ), 


4  V55  az^ 

where  i/’is  the  interface  value  of  m 


(2.30a) 


Similarly,  for  the  matrix  material,  if  we  multiply  Eq.  (2.9)  by  (  r  -r^  \  substitute  from 
Eq.  (2.21b)  and  integrate  by  parts  according  to  Eq.  (2.1  lb)  we  get 


(2-30b) 

4  AW55  5z 

where  we  have  used  the  continuity  relation  required  by  Eq.  (2.10d),  namely 


w*(z, t)  =  u^{z,ry,t)^  u„ (z, r, , /) . 


(2.31) 


Here  Q  =  +  3  +  21n  Wy),  is  a  function  of  the  volume  fractions  only.  By 


eliminating  u*  from  Eqs.  (2.30a)  and  (2.30b)  we  get  the  important  relation: 


-  -  . 


(2.32a) 
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where 


^,=i(I-0  and  =  2-).  (2,32b) 

4  4/55  m,. 

Next,  we  turn  our  attention  to  the  radial  and  tangential  stresses.  Subtracting  Eq. 
(2.8)  from  Eq.  (2.7)  we  get 

(2,33) 

or  r 

Specializing  this  equation  to  the  fiber  using  the  approximation  of  Eq.  (2.20a)  yields 

=Ggf,  (2.34a) 

as  one  expects  in  the  static  limit.  Specializing  the  same  equation  to  the  matrix,  using  Eq. 
(2.21a),  gives 


^rm  ^9m 


(2.34b) 


Following  the  above  procedure,  if  we  substitute  from  Eqs.  (2.20b)  and  (2.34a) 

into  Eq.  (2.2),  multiply  the  resulting  equation  by  r'^  and  integrate  according  to  Eq.  (2. 11a), 
we  get 


♦  _  r^  .dA  ^  n 

a.  -cr_r+— (— — Of — 5-)  =  0. 

4  ^az  de 


(2.35a) 


Similarly,  substituting  from  Eqs.  (2.21b)  and  (2.34b)  into  Eq.  (2.2),  multiplying  the 


resulting  equation  by  •”  r^  )  and  integrating  as  per  Eq.  (2. 1  lb)  we  obtain 


♦  _  2m...  ,  nQ.dA  dB 


(2.35b) 


where 
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a-/ (2,0  =  CT^fiz,r„t)  =  cr^„{z,r„t),  (2.36) 

is  the  common  interfacial  radial  stress. 

Eliminating  <j*  fromEqs.  (2.35a)  and  (2.35b)  yields 

_  _  r^Qr^A  d'^B^  t  ,  i  ^D 

a,r-(Jr„  =^-(— -P/- — - Pm — t) - ^(^m  +  « /•  m  77  ^  )5 .  (2.37) 

tn 

So  far  all  required  interface  conditions  have  been  identically  satisfied.  We  now 
need  to  eliminate  and  ^  from  the  above  equations.  In  order  to  do  this,  we  specialize 

the  constitutive  equation  (2.7)  for  the  fiber  and  matrix  by  substituting  into  it  from  Eqs. 
(2.20a)  and  (2.21a),  respectively  to  get 


^rf  =  (/22  +  fl^B  +  /j2 


du. 


dz  ’ 


(2.38a) 


O'rm  = - -(^22+  ^23  +  777,2 


77. 


dz 


(2.38b) 


Averaging  Eqs.  (2.38)  according  to  Eq.  (2.11)  yields,  for  the  fiber  and  the  matrix, 
respectively 


<^rf  -  if  12  +  fl3)B  +  /i2 


du, 


dz 


(2.39a) 


—  f  '2'Tn^ ,  .  „  du„ 

<^rn.  = - ^(^”22  +^23 - ^In 77^ ) 5  +  777,2 

77„  77„  dz 


(2.39b) 


Eliminating  cr,.^  and  between  Eqs.  (2.37)  and  (2.39a,  b)  we  finally  get  the  following 
important  equation 


dz  dz’ 


(2.40a) 


where 
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#1= - ["m(/22  +/23)  +  "/('”22+'”23)  +  2/W44]  End  ^2  =^iPf  -QPm)  ■ 

m 

(2.40b) 

The  behavior  of  the  composite  system  is  completely  described  by  Eqs.  (2.28), 
(2.29),  (2.32)  and  (2.40).  As  seen,  this  system  consists  of  four  partial  differential 
equations  that  couple  the  two  longitudinal  averaged  displacements,  interfacial  shear  stress 
and  interfacial  radial  displacement.  In  the  next  chapter  these  equations  will  be  utilized 
and  manipulated  to  facilitate  the  analysis  of  several  important  problems  involving  both 
static  and  dynamic  loading,  in  the  presence  or  absence  of  three  types  of  damage;  fiber 
breaks,  matrix  cracks  and  fiber-matrix  debonding. 
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CHAPTERS 


APPLICATIONS  OF  THE  STRAIGHT  FIBROUS  MODEL 

3.1  Introduction 

As  mentioned  in  chapter  1,  it  is  the  theme  in  this  work  to  use  straight  fibrous 
composites  analysis  as  the  first  building  block  unit  for  the  more  complex  analysis  of 
undulated  fibrous  composites.  It  is  important,  therefore,  to  have  a  deeper  understanding  of 
the  active  stress  transfer  mechanisms  that  takes  place  upon  loading  this  class  of 
composites.  These  mechanisms  change  according  to  the  type  of  loading  and  whether  the 
composite  has  some  sort  of  internal  damage  in  one  of  its  constituents.  The  systems  of 
equations  derived  in  the  previous  chapter  for  both  the  laterally  constrained  and  laterally 
free  cases  can  be  successfully  used  to  study  these  different  situations  involving  straight 
fibrous  composites.  The  simplest  of  which  is  the  problem  involving  only  static  loading 
with  or  without  damage.  They  can  also  be  used  to  solve  the  more  general  problem  of 
harmonic  dynamic  loading  of  the  straight  fibrous  composite. 

3.2  Cases  of  Static  Loading 

The  system  of  equations  (2.28),  (2.29),  (2.32)  and  (2.40)  for  the  case  of  laterally 
constrained  composite  segment  derived  in  the  previous  chapter  can  be  applied  to  cases  of 
static  loading,  by  simply  setting  time-dependent  terms  equal  to  zero.  Upon  doing  so,  this 
system  of  differential  equations  can  be  simplified  and  combined,  after  some  algebraic 
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manipulations,  into  a  single  fourth  order  differential  equation.  This  equation  takes  the 
form: 


a,a^'+  +  04  =  0  (3.1) 

where  the  prime  indicates  didz  and  a,-,  i=l,  2,  3,  are  functions  of  the  fiber  and  matrix 
material  properties  and  volume  fractions,  defined  as: 


'  2D 

1—^2  /„  ^12  I  /l2  \ 

T~Tn^  /T" 

2  ^11  fn 


^3 


n.  O  '  /„ 

with  ^4  depending  on  the  external  axial  loading  as 


P  „  2m, 2  .  m,2  M 

a. - P  + 

m 


(3.2a) 

(3.2b) 

(3.2c) 


11 


f\\ 


(3.2d) 


In  the  above  definitions,  D  is  given  by 

1  m  ^  f  ^ 

D  = — [n^  (m22  +  m23 )  +  n„  (/aa  +  /23 )  +  2^44  "  2(« /  — ^  +  n„  -f-)] . 

"111  /ii 

It  is  our  intent  to  drop  the  overscore  on  the  axial  stresses  and  displacements  from  now  on, 
since  no  more  averaging  is  going  to  take  place  in  the  following  chapters. 


3.2.1  Laterally  Free  System 

Laterally  free  systems  are  important  since  they  represent  a  widely  common  class 
of  composites,  namely  thin  fibrous  composites.  Due  to  the  similarity  in  the  field 
equations,  associated  interface  conditions  and  most  of  the  symmetry  conditions  between 
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the  laterally  free  and  the  laterally  constrained  systems,  we  shall  not  attempt  to  repeat  the 
procedure  described  above  for  this  case.  It  is  sufficient  to  mention  that  for  the  laterally 
free  situation,  the  symmetry  condition  in  Eq.  (2.10  b)  is  no  longer  valid  and  should  be 
replaced  with  the  appropriate  stress  free  boundary  condition,  namely 
^rmU>^2>0  =  0.  (3.3) 

With  this  restriction,  the  analysis  identically  parallels  that  of  the  laterally  constrained 
case,  and  the  resulting  system  of  equations  can  be  reduced  into  a  single  fourth  order 
differential  equation  identical  to  equation  (3.1).  In  this  case,  however,  the  coefficients  a,- 
take  different  forms. 

The  fourth  order  differential  equation  (3.1),  together  with  the  set  of  coefficients 
(3.2)  or  their  alternative  form  for  the  laterally  free  case,  are  general  and  can  be  applied  to 
a  variety  of  situations.  These  include  situations  pertaining  to  the  loading  conditions  and 
also  to  the  presence  or  absence  of  damage  in  the  system.  Also  these  equations  can  be 
applied  to  finite  as  well  as  to  infinitely  extended,  along  the  fiber  direction,  composites. 
Loading  conditions  consisting  of  uniform  heating,  mechanical  loading  or  combinations  of 
both  are  possible.  Damage  in  the  system  includes  the  possibility  of  matrix  cracking  or 
fiber  break.  A  complete  study  of  these  situations  was  published  in  a  recent  paper  [24]. 

3.2.2  Damage  Free  case 

For  the  case  of  perfectly  bonded  damage-free  fibrous  composite  under 
axisymmetric  loading,  the  solution  for  the  fiber  normal  stress  is  obtained  by  solving 
equation  (3.1)  subject  to  the  boundary  conditions 

o,^(0)  =  Po/  ,o^{L)  =  P,^  ,o;/(0)  =  0  ,o;^(L)  =  0  (3.4) 
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where  Po/  and  are  the  fiber  end  stresses  calculated  from  the  composite  total  stress  and 
the  fiber  volume  fraction.  The  appropriate  equations  of  the  analysis  of  chapter  2  can  then 
be  used  to  obtain  the  distributions  of  the  remaining  stress  components,  namely  the  matrix 
normal  stress,  the  shear  stress  and  the  radial  normal  stress. 

This  procedure  is  applied  for  the  case  shown  in  figure  3.1.  The  length  of  the 
composite  segment  is  chosen  as  L  =  8ri  and  the  fiber  has  a  volume  fraction  of  0.444.  The 
stresses  at  both  ends  are  taken  to  be  the  same  in  the  isotropic  fiber  and  matrix.  The  silica 
fiber  and  the  aluminum  alloy  matrix  have  the  following  values  for  Young’s  modulus  and 
Poisson’s  ratio: 


_ L _ ^ 

Figure  3.1:  Geometry  and  loading  of  the  case  of  a  finite  length  undamaged  composite. 
Pj,=73GPa  ,  =0.177  , 

P  =70.5GPa  ,  v„=  0.345. 

Figures  3.2,  3.3  and  3.4  show  the  resulting  distributions  of  the  fiber  normal  stress, 
interfacial  shear  stress  and  interfacial  radial  stress,  normalized  with  respect  to  the 
composite  total  stress  P.  Symmetry  of  the  solutions  are  evident  in  both  cases  of  laterally 


-  24 


free  and  laterally  constrained  segments.  In  the  absence  of  damage,  a  general  symmetric 
loading  along  the  fiber  direction,  namely  at  z  =  0  and  z  =  L  will  insure,  as  a  minimum,  the 
vanishing  of  shear  stress  and  normal  displacements  and  at  the  center  of  the  system 

i.e.,  at  z  =  U2.  The  effect  of  the  lateral  constraint  is  seen  to  increase  the  value  of  the 
interfacial  radial  stress. 


Figure  3.2:  Comparison  of  the  fiber  normal  stress  distribution  in  both  laterally  free  and 
laterally  constrained  systems  for  the  case  of  finite  undamaged  fiber-matrix  composite 
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z/r^ 


Figure  3.3:  Comparison  of  the  interfacial  shear  stress  distribution  in  both  laterally  free 
and  laterally  constrained  systems  for  the  case  of  finite  undamaged  fiber-matrix  composite 


Figure  3.4:  Comparison  of  the  interfacial  radial  stress  distribution  in  both  laterally  free 
and  laterally  constrained  systems  for  the  case  of  finite  undamaged  fiber-matrix  composite 
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It  is  important  to  emphasize  the  fact  that  despite  the  integration  operations  during 
the  analysis  in  chapter  2,  we  can  still  recover  the  distribution  of  stress  components  as 
functions  of  z  and  r.  For  example  the  shear  stress  distribution  in  both  the  fiber  and  matrix 
for  the  present  undamaged  fiber-matrix  composite  is  shown  in  figure  3.5.  Positive  and 
negative  values  of  the  stress  are  given  in  shades  of  white  and  black,  respectively. 


0  0.25  0.5  0.75  1  1.25  1.5 

z  (cm) 

Figure  3.5:  Distribution  of  the  shear  stress  in  an  undamaged 
fiber-matrix  composite  under  axial  loading. 

3.2.3  Broken  fiber  case 

For  cases  involving  a  fiber  break,  say  at  z  =  0  the  shear  stresses  and  and 
the  matrix  normal  displacement  vanish  at  that  location.  The  applicable  boundary 


conditions  are: 


In  all  the  numerical  illustrations  involving  damaged  composites,  arrangements 
studied  by  McCartney  [10]  have  been  used.  However,  these  arrangements  have  been 
studied  for  both  types  of  lateral  constraint.  Curves  belonging  to  the  laterally  free  system 
can  be  directly  compared  with  those  presented  in  McCartney’s  paper. 

As  a  numerical  example,  we  treat  the  situation  of  a  broken  fiber  in  a  perfectly 
bonded  fiber-matrix  infinite  system  in  the  absence  of  thermal  loading  as  shown  in  fig.  3.6. 
This  system,  which  was  described  in  [10],  section  8.1,  is  assumed  to  be  infinitely  long,  so 
that  the  fiber  end  stress  Pi/  is  equal  to  its  remote  value  based  on  an  equi-strain 
assumption,  namely  .  The  silica  fiber  and  the  aluminum  alloy  matrix  have  the  same 

properties  of  the  illustration  in  the  previous  sub-section  3.2.2.  For  the  arbitrarily  selected 

value  of  strain,  £  =  01 ,  the  remote  normal  stresses  have  the  following  values 

cr^^”  =7.3484  GPa  ,  =  7.0969  GPa  for  the  laterally  free  system,  and 

=8.333 IGPa  ,  =  12.260  GPa  for  the  laterally  constrained  system. 


Figure  3.6:  Geometry  and  loading  of  the  case  of  broken  fiber. 
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Solutions  for  this  case  are  obtained  from  Eqs.  (3.1)  with  coefficients  obtained 
from  (3.2),  subject  to  the  boundary  conditions  (3.5).  Results  of  the  fiber  normal  stress 
and  the  matrix  normal  stress  zie  plotted  and  compared  in  Fig.  3.7  for  both  the 

laterally  free  system  and  the  laterally  constrained  system.  In  the  absence  of  the  break,  the 
normal  stresses  in  both  the  fiber  and  matrix  are  constant  and  equal  to  their  asymptotic 

values  <7^“  and  aj°  as  predicted  by  the  Lame’  solutions.  In  the  presence  of  the  fiber 

break,  however,  this  figure  clearly  shows  the  normal  stress  transfer  between  the  fiber  and 
the  matrix.  From  this  figure  it  is  evident  that  the  matrix-fiber  load  transfer  takes  place  in 
the  vicinity  of  the  fiber  break  and  is  almost  complete  within  a  distance  of  2rj . 

The  interfacial  shear  stress  o*  and  the  interfacial  radial  stress  or*  for  the  above 

case  are  normalized  with  respect  to  the  remote  normal  matrix  stress  a ^  and  are  plotted 

and  compared  in  Figs.  3.8  and  3.9  for  both  lateral  constraint  cases.  Once  again,  both 
figures  clearly  show  that  most  of  the  activities  take  place  in  the  vicinity  of  the  break  and 
that  the  respective  stresses  asymptotically  approach  those  pertaining  to  the  damage  free 
systems  as  predicted  by  the  Lame’  solutions.  Figure  3.8  also  shows  that  the  condition  of 
zero  interfacial  shear  stress  at  the  fiber  breakage  location  is  exactly  satisfied  and  that  the 
location  of  the  maximum  value  of  a*  is  closer  to  the  location  of  the  broken  fiber  surface 
in  the  laterally  constrained  case  than  that  in  the  laterally  free  case.  The  maximum  value  of 
o*  itself  is  considerably  higher  in  the  former,  a*  is  seen  in  Fig.  3.9  to  be  dramatically 
affected  by  the  lateral  constraint.  This  is  perhaps  due  to  the  fact  that  the  laterally 
constrained  system,  in  which  the  lateral  strain  vanishes,  creates  large  radial  stresses  at  the 
outer  boundary  and  at  inner  locations  in  the  system. 
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Figure  3.7:  Comparison  of  the  normal  stress  distributions  in  the  fiber  and  the  matrix  for 
both  laterally  free  and  laterally  constrained  systems  in  the  presence  of  broken  fiber  at  z=0. 


Figure  3.8:  Comparison  of  the  interfacial  shear  stress  in  both  laterally  free  and  laterally 
constrained  systems  in  the  presence  of  broken  fiber  at  z  =  0. 
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Figure  3.9:  Comparison  of  the  interfacial  radial  stress  in  both  laterally  free  and  laterally 
constrained  systems  in  the  presence  of  broken  fiber  at  z  =  0. 

Finally,  the  effect  of  a  single  mid-span  fiber  break  on  the  shear  stress 
redistribution  in  the  initially  undamaged  composite  of  section  3.2.2  is  shown  in  figure 
3.10.  The  complete  antisymmetry  indicates  that  the  condition  of  zero  shear  stress  at  the 
break  surface  is  completely  satisfied. 

3.2.4  Cracked  matrix  case 

Opposite  to  the  case  of  the  broken  fiber,  cases  involving  matrix  cracking 
transverse  to  the  fiber  direction,  say  at  z  =  0,  will  exhibit  the  vanishing  of  the  fiber  normal 
displacement  at  that  location.  Here  the  applicable  boundary  conditions  then  become 
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Figure  3.10:  Distribution  of  the  shear  stress  in  a  fiber-matrix 
composite  under  axial  loading  with  a  single  fiber  break. 

o^(0)  =  Po/n^  ,a,^(L)  =  P^  ,o;^(0)  =  0  ,o'^{L)  =  Q  (3.6) 

As  an  illustration  we  treat  the  problem  of  a  cracked  matrix  in  a  perfectly  bonded 
fiber-matrix  infinite  system  with  thermal  loading,  as  shown  in  figure  3.11.  Again,  this 
system,  which  was  described  in  [10],  section  8.2,  is  assumed  to  be  infinitely  long.  The 
fiber  volume  fraction  is  0.4,  and  the  applied  strain  is  0.8%.  The  overall  system  is 
subjected  to  a  temperature  T  =  -550®C  below  the  stress  free  temperature.  Thermal  loading 
provides  no  more  difficulty  in  the  analysis,  since  its  effect  follows  that  of  the  external 
mechanical  loading  P,  and  thus  adds  to  its  term  in  equation  3.1.  The  fiber  and  matrix 
material  properties  are  isotropic  and  have  the  following  values  for  Young’s  modulus, 
Poisson’s  ratio  and  the  coefficient  of  thermal  expansion 
jE:^=220  GPa,  =0.1  ,  =5.5x10"^ /°C 

£„=70  GPa  ,  v„=0.2  ,  flr„=10-'/°C 
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Figure  3.1 1:  Geometry  and  loading  of  the  case  of  a  cracked  matrix  with  heating. 

Solutions  for  this  case  are  obtained  from  Eq.  (3.1)  subject  to  the  boundary  conditions 
(3.2).  The  remote  normal  stress  values  for  this  case  are  calculated  as 

=  1 .5866  GPa,  0.678 1  GPa  for  the  laterally  free  system,  and 

(jjf  =  1.7361  GPa,  =  0.9498  GPa  for  the  laterally  constrained  system. 

The  fiber  normal  stress,  the  interfacial  shear  stress  and  the  interfacial  radial  stress  are 
presented  and  compared  in  Fig.  3.12,  3.13  and  3.14.  These  figures  exhibit,  to  a  large 
degree,  the  overall  features  of  the  corresponding  figures  3.7,  3.8  and  3.9  in  both  regions 
of  maximum  activities  and  asymptotic  values.  It  is  worth  noting  that  the  lateral  constraint 
changed  the  nature  of  the  distribution  of  the  interfacial  radial  stress.  Instead  of  being 
totally  compressive  except  near  the  matrix  crack,  it  is  tensile  all  over  in  the  case  of 
laterally  constrained  system. 
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Figure  3.12:  Comparison  of  the  fiber  and  matrix  normal  stress  in  both  laterally  free  and 
laterally  constrained  systems  for  the  case  of  a  cracked  matrix  at  z  =  0  with  e  -  0.8%. 


Figure  3.13:  Comparison  of  the  interfacial  shear  stress  in  both  laterally  free  and  laterally 
constrained  systems  for  the  case  of  a  cracked  matrix  at  z  =  0  with  e  =  0.8%. 
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Figure  3.14:  Comparison  of  the  interfacial  radial  stress  in  both  laterally  free  Mid  laterally 
constrained  systems  for  the  case  of  a  cracked  matrix  at  z  =  0  with  £  =  0.8%. 

3.3  Cases  of  Dynamic  Loading 

We  already  established  in  the  previous  sections  of  this  chapter  that  the  analyses  of 
the  laterally  free  case  and  the  laterally  constrained  case  are  parallel.  We  will  choose  only 
one  of  them  for  any  single  application,  with  the  understanding  that  similar  argument  can 
be  established  for  the  other.  As  far  as  dynamic  loading  is  concerned,  we  will  deal  only 
with  the  laterally  constrained  case,  discussed  in  Chapter  2.  The  behavior  of  a  laterally 
constrained  composite  system  is  therefore  described  by  the  four  coupled  equations  (2.28), 
(2.29),  (2.32)  and  (2.40).  Unlike  the  case  of  static  loading,  these  equations,  in  their  full 
form  can  not  be  reduced  to  a  single  equation  similar  to  equation  (3.1),  and  have  to  be 
dealt  with  as  is.  In  order  to  do  so,  we  assume  a  steady  state  harmonic  solution  for  the 
system  as 


■  35 


Uj  =  Xjg 

(3.7a) 

■V  ^CR+ita 
=  ^2^ 

(3.7b) 

(3.7c) 

B=X,e 


OZ+iOM 


(3.7d) 


where  X,. ,  i  =1..4  are  amplitudes  yet  to  be  determined  as  part  of  the  solution,  cd  is  the 
circular  frequency  and  oris  an  unknown  parameter. 

Substituting  these  expressions  into  the  system  of  equations  (2.28,  2.29,  2.32,  2.40)  and 


rearranging  in  matrix  form,  we  get  for  the  case  of  isotropic  fiber 


rif  (Aj-  +  2tij  )a^  +  TijPjQ)^ 

0 

Irij 

2nfA^a 

X/ 

0 

”m  (■^m  ^mPm^ 

—  2n^ 

-2n^A^a 

X^ 

-1 

1 

-Vi 

r],a 

-Xja 

T},a 

.^4. 

(3.8) 


For  nontrivial  solutions,  the  determinant  of  the  left  hand  side  must  be  zero.  Upon  doing 
this  we  get  six  values  for  the  exponent  a .  We  have  a  complete  set  of  coefficients  X, , 
1=1.. 4  corresponding  to  each  value  of  o  .  At  this  point  of  the  solution  only  their  ratios  can 
be  determined  and  therefore,  we  chose  to  solve  for  X2,  X3  and  X4  in  terms  of  Xj. 

Using  the  method  of  superposition,  the  formal  solutions  for  the  displacement 


,  A,S)=  2^,  ,V„ ,  V„  >f , 

S  (a.  .  ,  D,>.,  ^  . 

9=1 

where 


(3.9) 

(3.10) 


Dig  ~  +  2/^/  )^g  ’ 


36 


The  external  loading,  as  well  as  the  two  damage  types  are  now  imposed  as  boundary 
conditions  to  the  problem. 

3.3.1  Undamaged  fiber  matrix  Composite 

In  the  absence  of  damage,  the  applicable  boundary  conditions  at  both  ends  of  the 
composite  can  be  summarized  as 

a^{-Ll2)  =  Pe‘‘“  ,  £r^(-L/2)  =  Pe"“  ,  <7*(-L/2)  =  0  (3.11a) 

(T^(LI2)  =  Pe“^  ,  (7^(112)==  Pe*‘“  ,a\LI2)  =  0  (3.11b) 

Where  fJo ,  ,  i  =  1,2  are  the  fiber  and  matrix  normal  stresses  at  both  ends.  Since  all 

load  amplitudes  are  chosen  equal,  the  deformation  of  the  system  is  symmetric  about  its 
center  z=0  and  an  equivalent  set  of  boundary  conditions  will  then  be 
M^(0)  =  M„(0)  =  (T*(0)  =  0  and  also  either  of  equations  (3.11a)  or  (3.11b)  is  to  be 
satisfied.  Solutions  for  this  case  are  obtained  from  Eqs.  (3.9)  and  (3.10),  where 

values  are  obtained  from  satisfying  the  boundary  conditions  (3.11). 

The  effect  of  the  loading  frequency  at,  in  krad/s,  on  the  distributions  of  the 
amplitudes  of  the  fiber  average  normal  stress  ,  the  matrix  average  normal  stress 

and  the  interfacial  shear  stress  o*  is  displayed  in  Figs.  3.15, 3.16  and  3.17,  respectively. 
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Figure  3.15:  Variation  of  the  amplitude  of  the  fiber  normal  stress  distribution  with  the 
loading  frequency  for  the  case  of  fiber  pull-out. 

In  these  figures  the  fiber  and  matrix  properties  are  the  same  as  those  used  in  the  cases  of 
static  loading,  and  only  the  fiber  is  loaded,  to  simulate  the  important  case  of  fiber  pull¬ 
out.  In  all  figures  we  notice  that  most  stress  variations  take  place  in  the  vicinity  of  the 
outer  edges  near  the  applied  loads.  For  each  frequency,  moving  from  the  loaded  end 
inward,  we  see  that  the  fiber  normal  stress  decreases  to  its  minimum  value  at  the  center 
and  that  of  the  matrix  increases  to  the  limiting  value  at  the  same  location.  The  interfacial 
shear  stress  starts  and  ends  with  zero  values  at  the  outer  and  center  locations  but  has  a 
sharp  variation  close  to  the  loaded  ends.  For  normal  stresses,  it  is  seen  that  an  increase  in 
frequency  leads  to  a  pronounced  stiffening  of  the  material  and  hence  to  an  increase  in  the 
values  of  these  stresses.  This  is  in  contrast  with  the  interfacial  shear  stress,  which  seem  to 
be  rather  insensitive  to  frequency  variations. 
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Figure  3.16:  Variation  of  the  amplitude  of  the  matrix  normal  stress  distribution  with  the 
loading  frequency  for  the  case  of  fiber  pull-out. 


The  above  observations  suggest,  at  least  for  the  fiber,  that  fiber  breaks,  if  occur, 
will  most  probably  take  place  at  the  center  of  the  slab.  For  this  reason,  in  the  next  two 
subsections  we  discuss  the  possibilities  of  a  matrix  crack  or  a  fiber  break  at  the  center  of 
the  slab  separately. 


3.3.2  Broken  Fiber  Case 

For  the  case  of  a  broken  fiber  at  z  =  0,  the  shear  stress  in  both  media,  the  normal 
stress  in  the  fiber  and  the  displacement  in  the  matrix  all  vanish  at  the  break  location. 
Here,  the  applicable  boundary  conditions  then  become 
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Figure  3.17:  Variation  of  the  amplitude  of  the  interfacial  shear  stress  distribution  with  the 
loading  frequency  for  the  case  of  fiber  pull-out. 


o^(0)  =  0  ,  M„(0)  =  0  ,  ct\0)  =  0  (3.12a) 

(7^iLI2)  =  Pe““  ,  dr^(L/2)  =  Pe"“  ,  cr*(L/2)  =  0  (3.12b) 

Under  these  boundary  conditions,  we  solve  for  the  appropriate  amplitudes 

from  Eqs.  (3.9)  and  (3.10).  The  effect  of  the  load  frequency  on  the  distributions  of  the 
amplitudes  of  the  fiber  average  normal  stress  ,  the  matrix  average  normal  stress 

and  the  interfacial  shear  stress  o*  is  displayed  in  Figs.  3.18, 3.19  and  3.20,  respectively. 
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Figure  3.18:  Variation  of  the  amplitude  of  the  average  fiber  normal  stress  distribution 
with  the  loading  frequency  for  the  case  of  a  broken  fiber  at  z  =  0. 

In  these  figures,  once  again,  the  results  are  presented  for  the  same  frequencies 
used  in  the  previous  situation.  All  figures  clearly  show  comparative  increase  in  the 
stresses  for  increasing  frequency.  All  stresses  are  seen  to  be  more  or  less  similar  to  those 
corresponding  to  the  undamaged  case  in  the  vicinity  of  the  edge  loads.  Pronounced 
changes  in  the  stress  pattern  take  place  in  the  vicinity  of  the  break.  This  is  due  to  the  large 
transfer  of  the  load  between  the  matrix  and  the  fiber.  Beside  such  a  transfer  between  the 
normal  stress  components,  a  pronounced  variation  of  the  interfacial  shear  stress  takes 
place  close  to  the  break.  It  should  be  noted  that  the  crack  opening  can  be  obtained  from 
Eq.  (2.32)  by  setting  =0  and  A=  0,  resulting  in 
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Figure  3.19:  Variation  of  the  amplitude  of  the  matrix  normal  stress  distribution  with  the 
loading  frequency  for  the  case  of  a  broken  fiber  at  z  =  0. 
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Figure  3.20:  Variation  of  the  amplitude  of  the  interfacial  shear  stress  distribution  with  the 
loading  frequency  for  the  case  of  a  broken  fiber  at  z  =  0. 
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3.3.3  Matrix  Cracking 

For  cases  involving  matrix  cracking  transverse  to  the  fiber  direction,  say  at  z  =  0, 
the  shear  stresses  and  normal  stress  O2  the  matrix  and  the  fiber  normal 

displacement  SJ  vanish  at  that  location.  Within  the  context  of  the  present  modeling,  the 
applicable  boundary  conditions  in  the  presence  of  the  matrix  crack  can  be  summarized  as 
Oz.(0)  =  0  ,  «/0)  =  0  ,  (T*(0)  =  0  (3.14a) 

(T^iL/2)  =  Pe'“  ,  <7^(L/2)  =  Pe'“  ,  <7*(L/2)  =  0  (3.14b) 

Solutions  for  this  case  are  obtained  from  Eqs.  (3.9)  and  (3.10),  in  a  manner 
identical  to  that  of  the  broken  fiber.  Therefore  no  further  discussion  of  this  case  in 
necessary.  As  a  byproduct  of  the  analysis,  the  matrix  crack  opening  can  be  easily 
calculated  from  equation  (2.32)  by  merely  setting  =  0  and  A  =  0,  resulting  in 

«2(0)  =  -^,^  (3.15) 

The  formulation  of  section  3.3,  as  well  as  its  application  to  several  situations,  is 
detailed  in  [26].  There  are  several  ways  to  compare  the  results  of  the  dynamic  loading 
formulation  with  experimental  data.  One  is  to  compare  the  phase  velocities  of 
propagating  waves  at  different  frequencies  (dispersion  curves).  This  is  done  in  [27]  for 
the  case  of  a  titanium  alloy  matrix  reinforced  with  three-layer  concentric  cylindrical  SiC 
fibers.  Each  fiber  is  about  140  /nn  in  total  diameter  and  consists  of  a  carbon  core  and  a 
SiC  shell  separated  from  the  matrix  by  another  3-/fln-thick  carbon  layer.  The  mechanical 
properties  of  all  the  constituents  are  presented  in  table  3.1.  The  effective  fiber  anisotropic 
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properties  are  calculated  from  these  individual  isotropic  properties  by  a  repetitive 
application  of  the  procedure  outlines  in  [1],  which  can  handle  single  material  fiber  and 
matrix.  The  frequency  dependence  of  the  fundamental  mode  is  plotted  in  figure  3.21  and 
compared  with  the  exact  solution  and  with  the  experimental  data  of  Huang  and  Rokhlin. 
As  can  be  seen,  significant  improvements  are  achieved  by  the  improved  mixture  model 
when  compared  with  the  experimental  data  as  well  as  with  the  exact  solution. 


Table  3.1:  Material  properties  and  dimensions  for  the  fiber  and  matrix  materials. 


Material 

Density 

Cii 

Cl2 

Ci3 

C44 

C55 

Radius 

kg/m^ 

GPa 

GPa 

GPa 

GPa 

GPa 

jjm 

Isotropic  carbon 

1700 

49 

17 

17 

16 

49 

18 

SiC 

3200 

446 

92 

92 

177 

446 

68 

Interphase  carbon 

2100 

31 

21.8 

21.8 

4.6 

31 

71 

Titanium 

5400 

193 

103 

103 

45 

193 

132.3 

In  this  chapter  we  have  shown  the  versatility  of  the  elasticity-based  formulation  of 
the  micromechanical  behavior  of  straight  fibrous  composites.  The  modeling  was  found  to 
be  applicable  to  a  wide  range  of  important  problems,  and  its  validity  was  established  by 
its  ability  to  reproduce  available  published  results  [10].  Its  superiority  over  these  previous 
modelings  lies  in  its  ability  to  solve  problems  that  could  not  be  solved  using  these  models 
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Figure  3.21:  Fundamental  mode  comparisons  of  the  suggested  mixture  model, 
the  exact  solution  and  experimental  data. 


in  particular  problems  involving  finite  length  representative  cells  and  transversely 
isotropic  constituents. 
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CHAPTER 4 


DEVELOPMENT  OF  THE  MODEL  FOR  UNDULATED  FIBROUS 
COMPOSITES 

4.1  Introduction 

The  model  developed  in  chapter  2  for  straight  fibrous  composites  will  be  extended  to 
cases  of  undulated  fiber  reinforcement  under  axial  loading.  This  takes  the  analysis  one 
step  further  towards  modeling  plain  woven  composites.  For  simplicity,  we  assume  the 
undulation  to  be  restricted  to  a  single  plane,  say  the  x-y  plane.  The  procedure  is  outlined 
as  follows;  First,  for  a  given  geometric  undulation,  we  identify  and  analytically  describe 
local  tangents  to  the  fiber.  Second,  we  transform  the  global  coordinates  and  applied  loads 
to  local  coordinate  systems  consisting  of  the  tangent  directions  and  their  normals  in  the 
plane  of  undulation.  Third,  we  adapt  the  results  obtained  for  straight  fibers  to  "straight" 
segmented  fibers  along  the  tangents.  The  local  stresses  obtained  vmder  this  procedure  are 
then  supplemented  with  the  stress  contributions  that  inherently  rise  in  any  oriented 
direction  with  respect  to  the  loading  direction.  The  resulting  local  stresses  are  found  to 
compare  favorable  with  finite  element  calculations. 

4.2  Geometric  Model 

A  typical  undulated  fiber  reinforced  composite  is  shown  in  figure  4.1.  The 
periodicity  of  the  medium  allows  us  to  model  it  as  a  concentric  cylindrical  system 
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subjected  at  its  outer  boundary  to  vanishing  shear  stress  and  radial  displacement  as 
illustrated  in  figure  4.2.  The  fiber  centerline  is  assumed  to  follow  the  periodic  function 


Figure  4. 1 :  Undulated  fibrous  composite  panel. 


fix)  =  A  cos 


^  27tn{x-  E)' 


,n  =  \,l,  3,. 


(4.1) 


where  is  the  panel  length.  At  any  global  location  x,  the  local  inclination  angle  0  is 
then  defined  as 


6{x)  =  tan 


-I 


dx  J 


(4.2) 


Along  the  tangent  direction  at  the  angle  0 ,  a  “local”  coordinate  z  is  chosen  which  relates 
geometrically  to  the  global  coordinate  x  as 


z  =  x/  cos^(x) 


(4.3) 


-47 


Figure  4.2:  Representative  unit  cell. 


4.3  Derivation  of  the  Undulated  Case  from  the  Straight  Fibrous  Case 

The  representative  unit  cell  of  figure  4.2  is  assumed  to  be  loaded  symmetrically  as 
shown  in  figure  4.3.  The  loads  consist  of  axial  as  well  as  transverse  stress  components. 
The  axial  boundary  stresses  are  arbitrarily  assigned  the  values  Px\  and  Pjc2,  for  the  fiber 
and  matrix,  respectively.  The  transverse  stress  Py  is  still  an  unknown  and  rises  as  a 
consequence  of  the  axially  loaded,  but  laterally  constrained,  systems.  Its  value  will  be 
approximated  and  calculated  based  upon  the  Lame'  type  solutions  of  the  system  (under 
axial  iso-strain  assumptions).  Based  on  the  iso-strain  analysis,  Py  takes  the  form 


P  = 


(4.4) 


Along  the  tangent  direction,  i.e.  at  the  angle  of  inclination  6,  the  fiber  and  matrix 
boundary  stresses  transform  to 


P,  =  cos^  9-\-Py  sin^  $ 


(4.5a) 


Pi  -  Pxi  0^-Py  sin  6 


(4.5b) 


Figure  4.3 ;  Characterization  of  loads. 


The  fundamental  assumption  is  that  the  local  stresses,  at  any  location  x,  constitute 
a  superposition  of  two  components.  With  reference  to  figure  4.4,  the  first  is  obtained 
from  adapting  the  solutions  of  the  straight  fiber  to  directions  along  the  tangents  and  the 
second  is  the  inherent  stresses  that  rise  in  the  matrix  in  directions  that  are  inclined  to  the 
loading  directions.  The  first  component  is  more  difficult  to  obtain  and  requires  the 
solution  of  the  field  equations  for  the  straight  fiber  case  in  cylindrical  coordinate  system. 
This,  however,  was  treated  in  chapter  2. 

The  equations  developed  in  chapter  2  can  be  applied  to  any  straight  fiber  system, 
having  any  length  and  subjected  to  prescribed  axial  stresses  at  its  ends.  In  particular  this 
equation  can  be  updated  and  used  for  all  segments  along  the  undulated  fiber  by  choosing, 
for  a  given  0fxj,  the  value  of  /*  as  constructed  in  Eq.(4.5),  and  the  system  length  of 
L  =  LJ  cos  0(x) .  This  constitutes  the  first  part  of  the  stress  distribution.  The  second  part 
is  the  inherent  stress  distribution  that  rises  on  planes  that  are  inclined  at  the  0-  direction 
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from  the  loading  direction.  Thus  the  final  interfacial  shear  stress  distribution  takes  the 
form 

=  -{P^-P^s\n 9 CO&0 sin <p  (4.6) 

where  (f>  is  the  azimuthal  angle  that  runs  counter-clock  wise  along  the  circumference  of 
the  fiber.  Thus,  it  defines  the  x-r  plane  that  is  inclined  to  the  plane  of  undulation.  We 
note  that  the  maximum  contribution  of  the  inherent  stress  components  exists  at  ^  =  ±90° . 
These  are  the  directions  of  the  plane  of  undulation.  <f>  =  (f  and  180°  define  the  directions 
that  are  normal  to  the  plane  of  undulation  where  the  inherent  stresses  vanish.  It  is 
obvious  that  the  undulated  fiber  results  reduce  to  the  corresponding  case  of  straight  fibers 
by  merely  setting  the  amplitude  of  undulation  to  0. 

In  the  following  chapter  we  provide  comparisons  of  the  results  of  this  formulation 
with  those  of  finite  element  for  a  variety  of  situations.  We  also  discuss  the  applicability 
of  the  formulation  to  cases  involving  different  types  of  damage  and  static  as  well  as 
dynamic  loading. 
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CHAPTERS 


APPLICATIONS  OF  THE  UNDULATED  FIBROUS  MODEL 

5.1  Application  to  Cases  of  Static  Loading 

The  solution  strategy  described  in  chapter  2  for  the  case  of  straight  fibrous 
composites  and  extended  in  chapter  4  for  undulated  fibrous  composites  can  be  applied  to 
cases  of  static  loading.  This  is  done  in  details  in  [28],  In  that  case,  the  four  coupled 
differential  equations  (2.28),  (2.29),  (2.32),  (2.40)  can  be  combined  into  a  single  equation 
identical  to  equation  (3.1).  The  longitudinal  coordinate  z,  however,  is  now  in  the  local 
tangent  direction  for  each  segment.  Loads  also  have  to  be  transformed  into  local 
coordinates  as  explained  in  chapter  4.  The  procedure  is  demonstrated  in  this  section  for 
two  test  cases  and  the  results  are  compared  to  those  of  the  finite  element.  The  test  cases 
consist  of  a  fiber-matrix  composite  with  eight  complete  fiber  sinosoidal  waves.  Material 
properties  are  the  same  as  those  used  in  [10]  and  in  the  examples  of  Chapter  3.  The 
undulated  fiber  has  a  length  to  amplitude  ratio  L^l  a  =  %  and  a  fiber  volume  fraction  of 
0.321. 

Finite  element  results  are  generated  using  ANSYS  package.  The  FE  code  source 
file  itself  is  generated  using  a  modified  MATHEMATICA  program  that  was  originally 
developed  by  Al-Huniti  [2].  This  program  requires  the  fiber  and  matrix  properties,  the 
length/fiber  ratio,  and  the  element  size.  Only  half  of  the  problem  is  solved  to  make  use  of 
the  symmetry  about  the  x-y  plane.  Both  tetrahedron  and  brick  elements  are  used  and  the 
problem  symmetry  about  the  _y-axis  is  utilized.  Overall,  278  higher  order  tetrahedron 
elements  were  used  to  model  the  fiber  and  the  matrix  volumes.  The  fiber-matrix 
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composite  displacements  are  constrained  in  the  lateral  direction  in  order  to  satisfy  the 
required  condition  V2(z,rj)  =  0 ,  except  for  the  surface  x  =  1 2,  on  which  the  external 

load  is  applied.  Figure  5.1  shows  the  finite  element  model  used  in  the  test  cases. 


In  the  first  test  case,  the  undamaged  fiber-matrix  composite  is  under  condition  of 
fiber  pull-out.  Therefore,  the  matrix  surface  at  both  ends  is  stress-  firee.  Results  of  the 
fiber  normal  stress  and  the  interfacial  shear  stress  a*  are  shown  in  Figs.  5.2  and  5.3.  In 
calculating  from  the  finite  element  results  at  specific  cross-sections  (results  denoted 
by  triangles  on  the  figure),  the  stress  tensor  component  in  the  fiber  direction  was 
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Figure  5.2:  Comparison  of  the  fiber  normal  stress  distribution  during  fiber  pull-out  in  an 
8  wave  fibrous  composite  calculated  by  the  proposed  model  and  finite  element  results. 


X  (cm) 

Figure  5.3:  Comparison  of  the  interfacial  shear  stress  distribution  during  fiber  pull-out  in 
an  8-wave  fibrous  composite  calculated  by  the  proposed  model  and  finite  element  results. 
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averaged  over  the  fiber  cross  sectional  area.  A  good  agreement  between  the  analytically 
predicted  values  and  the  finite  element  solution  can  be  seen.  As  expected,  the  effect  of 
undulation  is  a  maximum  at  locations  of  maximum  slope  0,  and  the  stress  distribution 
approaches  that  of  the  straight  fibrous  case  at  locations  of  zero  slope. 

In  the  other  test  case,  damage  is  introduced  in  the  form  of  a  fiber  break  at  z  =  0. 
The  composite  is  subjected  to  an  overall  strain  of  0.1  %  at  both  ends.  As  seen  from 
figures  5.4,  there  exists  a  very  good  agreement  between  the  predicted  distribution  for  the 


X  (cm) 

Figure  5.4:  Comparison  of  the  fiber  normal  stress  distribution  in 
the  presence  of  a  fiber  break  in  an  8  wave  fibrous  composite  under  axial 
loading  calculated  using  the  proposed  model  and  finite  element  results. 

fiber  normal  stress  and  the  finite  element  results.  Results  are  not  so  good  in  the  case  of 
interfacial  shear  stress,  as  seen  in  figure  5.5.  The  reason  is  that  the  calculation  of  a 
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involves  an  additional  differentiation  in  both  the  analytical  and  the  finite  element 
solution.  This  usually  gives  rise  to  some  deviation  in  the  results. 


Figure  5.5:  Comparison  of  the  interfacial  shear  stress  distribution  in 
the  presence  of  a  fiber  break  in  an  8  wave  fibrous  composite  under  axial  loading 
calculated  using  the  proposed  model  and  finite  element  results. 

5.2  Cases  of  Dynamic  Loading 

As  explained  in  chapter  4,  the  solution  method  discussed  in  chapter  2  can  be 
adopted  to  cases  of  undulated  reinforcement  if  the  local  coordinate  as  well  as  the  local 
transformed  loads  are  used.  As  illustrations  of  this  case,  we  once  more  revisit  the  two  test 
cases  of  the  previous  section  5.1.  Here,  however,  we  are  to  investigate  the  effect  of  a 
harmonic  type  of  loading  on  the  stress  distributions.  Figures  5.6,  5.7  and  5.8  show  that 
effect  for  the  case  of  fiber  pull-out  for  frequencies  below  the  resonance  frequency.  Both 
the  fiber  and  matrix  normal  stresses  increase  considerably  with  the  loading  frequency. 
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while  the  interfacial  shear  stress  distribution  is  almost  independent  of  it.  Figure  5.8  shows 
therefore,  that  the  progression  of  delamination,  which  is  caused  mainly  by  the  interfacial 
shear  stress  a ,  does  not  depend  on  the  loading  frequency  for  the  case  of  harmonic 
loading. 


X  (cm) 

Figure  5.6:  Variation  of  the  amplitude  of  the  fiber  normal  stress  distribution  with  the 
loading  frequency  for  the  case  of  undulated  fiber  during  fiber  pull-out. 

Next  we  consider  the  other  test  case  of  an  undulated  fibrous  composite  that 
include  broken  fibers  under  applied  harmonic  strain  of  amplitude  0.1%.  Figures  5.9,  5.10 
and  5. 1 1  show  how  the  amplitudes  of  the  various  stress  components  vary  with  the  loading 
angular  frequency.  Figures  5.9  and  5.10  show  that  the  normal  stresses  exhibit  a  trend 
similar  to  that  of  the  corresponding  figures  of  the  fiber  pull-out  case.  The  interfacial  shear 
stress  near  the  damage,  however,  depends  considerably  on  the  loading  frequency,  as  seen 
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Figure  5.7:  Variation  of  the  amplitude  of  the  matrix  normal  stress  distribution  with  the 
loading  frequency  for  the  case  of  undulated  fiber  during  fiber  pull-out. 


X  (cm) 

Figure  5.8:  Variation  of  the  amplitude  of  the  interfacial  shear  stress  distribution  with  the 
loading  frequency  for  the  case  of  undulated  fiber  during  fiber  pull-out. 
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Figure  5.9:  Variation  of  the  amplitude  of  the  fiber  normal  stress  distribution  with  the 
loading  frequency  for  the  case  of  undulated  fiber  with  a  break. 
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Figure  5.10:  Variation  of  the  amplitude  of  the  matrix  normal  stress  distribution  with  the 
loading  frequency  for  the  case  of  undulated  fiber  with  a  break. 
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in  figure  5.11.  Since  delamination  usually  starts  at  the  free  surface  of  a  damaged  fiber  or 
matrix,  this  suggests  the  dependence  of  the  delamination  propagation  rate  in  this  case  on 
the  loading  frequency. 

Another  important  variable  that  can  be  predicted  by  the  model  is  the  crack 
opening  displacement.  The  amplitude  of  the  opening  A  =  m/0)  is  calculated  as  per 
equation  (3.13)  and  normalized  with  respect  to  its  value  at  zero  loading  frequency  (static 
loading).  The  variation  of  this  normalized  displacement  with  the  loading  frequency,  up  to 
1.5  ftbs,  is  shown  in  figure  5.12.  Results  of  other  cases  of  the  behavior  of  undulated 
fibrous  composites  subject  to  harmonic  loading  are  published  in  [29]. 


X  (cm) 

Figure  5.1 1 :  Variation  of  the  amplitude  of  the  interfacial  shear  stress  distribution  with 
the  loading  frequency  for  the  case  of  undulated  fiber  with  a  break. 
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Figure  5.12:  Variation  of  the  amplitude  of  the  fiber  break  opening  displacement  with  the 
loading  frequency  for  the  case  of  undulated  fiber. 

With  these  two  test  cases,  the  analysis  of  the  basic  element  of  woven  composites, 
namely  the  undulated  fibrous  composite  is  concluded.  The  micromechanical  model 
provided  accurate  predictions  for  the  stresses  in  both  of  the  composite  constituents,  which 
are  responsible  for  constituent  failure,  as  well  as  for  the  interfacial  shear  stress, 
responsible  for  the  delamination  between  the  fiber  and  the  matrix.  These  predictions  will 
be  used  in  the  next  chapter  to  predict  the  progression  of  transverse  cracks  as  well  as  the 
extension  of  regions  of  delamination.  Degradation  of  the  composite  properties  that  takes 
place  will  accordingly  be  quantified. 
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CHAPTER  6 


DAMAGE  ANALYSIS  IN  FIBROUS  COMPOSITES 

6.1  Introduction 

The  term  damage  mechanics  is  usually  used  to  describe  the  study  of  the 
accumulation  of  damage  from  its  initiation  up  to  and  including  rupture.  Damage  is 
usually  an  ill-defined  term.  By  damage  here,  we  mean  a  local  micromechanical  failure  at 
the  fiber/matrix  level.  This  failure  occurs  through  one  of  several  mechanisms  including 
fiber  breaking,  matrix  cracking,  fiber-matrix  debonding  and  fiber  pull-out.  The  existence 
of  a  distribution  of  fiber  breaks  or  matrix  transverse  cracks  significantly  degrades  the 
stiffness  and  residual  strength  of  the  composite.  It  is  the  purpose  of  this  chapter  to  study 
such  an  effect  on  straight  fibrous  composites. 

In  this  chapter  we  will  use  our  basic  model  to  predict:  1.  whether  failure  takes 
place  in  an  undamaged  composite  panel  under  certain  static  loading,  and  2.  how  much  the 
successive  failure  affects  the  composite  stiffness  as  expressed  by  Young’s  modulus.  The 
following  two  sections  explain  the  procedure  to  achieve  these  goals.  Our  modeling  will 
adopt  the  maximum  stress  criterion  for  the  onset  of  failure. 

6.2  Damage  Accumulation  Prediction 

Consider  an  undamaged  straight  fibrous  composite  as  shown  in  figure  2.1.  Upon 
application  of  external  loading  P,  equation  (3.1)  can  be  solved  for  the  fiber  normal  stress 
distribution  in  the  representative  unit  cell.  From  the  equilibrium  of  the  internal  stresses  at 
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each  section  with  the  external  applied  load,  the  normal  stress  distribution  in  the  matrix  is 
calculated  from 


=  (61) 

Using  the  maximum  stress  criterion,  the  matrix  will  fail  when 

jr,-Afe[<7„)£0  (6.2) 

where  Xm  is  the  strength  of  the  matrix  material.  The  same  is  true  for  the  fiber  material. 
When  failure  takes  place  in  either  the  fiber  or  the  matrix,  the  unit  cell  length  should  be 
adjusted  accordingly  and  the  process  is  repeated  for  increasing  load  values. 


6.3  Evaluation  of  Young’s  Modulus 

The  procedure  to  calculating  the  longitudinal  stiffness  at  a  general  loading  step 
point  is  described  in  this  section  for  the  case  of  an  existing  fiber  break  in  the  system.  A 
parallel  argument  can  be  made  for  the  case  of  an  existing  matrix  crack.  As  mentioned  in 
the  previous  section,  the  first  generation  of  results  in  this  modeling  is  always  the  fiber 
normal  stress  as  obtained  from  solution  of  (3.1).  Next,  equation  (6.1)  is  used  to  determine 

the  matrix  normal  stress  distribution  .  Once  the  normal  stresses  in  the  fiber  and 
matrix  are  known,  the  integrated  constitutive  relation  (2.27)  can  be  used  to  evaluate  the 
distribution  matrix  longitudinal  strain,  resulting  in 


dz  m,, 


(6.3) 


Integrating  this  strain  over  the  length  of  the  unit  cell,  namely  fi'om  z  =  -Z,;/2  to  z  =  Z,//2 
provides  the  distribution  of  longitudinal  displacements  in  the  matrix  media  u^{z) .  The 
total  composite  Young’s  modulus  is  therefore  calculated  from 


(6.4) 


F  = _ - 

•4-  2»,(£,/2) 

The  procedures  provided  in  sections  6.2  and  6.3  are  better  illustrated  via  a  flow  chart. 
This  is  shown  in  figure  6.1 

Figure  6.2  demonstrates  the  process  of  damage  propagation,  in  terms  of  transverse 
crack  accumulation,  in  a  straight  fibrous  composite.  The  external  loading  is  normalized 
such  that  failure  occurs  when  the  loading  has  the  nominal  value  of  10  MPa.  The  fiber  and 
matrix  are  E-glass  and  Epoxy,  respectively.  The  mechanical  properties  of  these  two 
materials  are  obtained  from  reference  [24]  as 
Ej-  =  69  GPa,  Vj.  =  0.22 ,  Xj-  =  3450  MPa, 

E„  =  4.6  GPa,  v„  =  0.36,  X„  =  58.6 MPa 

From  the  figure,  it  is  clear  that  in  the  case  of  straight  fibrous  composites  the 

representative  unit  cell  is  divided  by  the  cracks  into  pieces  of  approximately  equal  length. 

This  is  because  the  normal  stress  distribution  between  two  cracks  has  a  maximum  value 
at  the  mid-point  between  these  cracks.  There  is  experimental  evidence  that  the  mean 
value  of  the  probabilistic  distribution  of  the  crack  location  is  close  to  the  middle  of  the 
segment.  Because  of  this,  we  are  able  to  define  the  crack  density  for  a  straight  fibrous 

composite  as  c  =  1/Z,.,  where  Z;  is  the  distance  between  two  adjacent  cracks  at  damage 
step  /. 
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Start  with  undamaged  unit  cell  of 
length  Li  subjected  to  a  small  load  Pi 


Figure  6. 1 :  Flow  chart  of  calculating  the  damage  accumulation  and  property  degradation. 
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Figure  6.2:  Damage  accumulation  in  an  E-glass/Epoxy  fibrous  composite. 
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The  degradation  curve  of  Young’s  modulus  of  this  composite  as  calculated  from 
equation  6.4  is  shown  in  figure  6.3  from  crack  initiation  (no  damage),  to  crack  saturation 
{Li  is  very  small).  The  value  of  Young’s  modulus  is  normalized  with  respect  to  its 
undamaged  value  calculated  fi’om  Lame  solution  as  73  GPa.  At  very  high  crack  density, 
the  value  of  Young’s  modulus  is  identical  to  the  case  without  the  fiber. 


Figure  6.3:  Degradation  curve  for  an  E-glass/Epoxy  Composite 
from  crack  initiation  to  crack  saturation. 

A  modeling  similar  to  that  presented  in  chapter  2  and  section  6.3  can  be  used  to 
study  the  problem  of  progressive  transverse  cracking  in  symmetric  laminated  composites. 
This  is  useful  as  most  of  the  experiments  are  conducted  on  the  commonly  used  laminated 
composites.  In  a  recent  publication  [30]  the  problem  of  transverse  cracking  in  the  90° 
layer  of  a  balanced  symmetric  [0°/90°]  laminate  was  studied.  The  E-glass/Epoxy  laminate 
was  loaded  in  the  0°  fibers’  direction.  A  comparison  of  the  predicted  variation  of  the 
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normalized  Young’s  modulus  with  the  experimental  measurements  showed  perfect 
agreement.  For  illustration,  this  comparison  is  reproduced  in  figure  6.4.  As  the  figure 
shows,  accuracy  of  the  predicted  degradation  is  maintained  all  the  way  from  crack 
initiation  (c=0)  to  the  crack  saturation  case. 


c  (  1/mm) 

Figure  6.4:  Comparison  between  the  predicted  distribution  of  the  residual 
Young’s  modulus  and  the  experimental  measurements. 

Having  established  the  procedure  of  predicting  the  stiffness  degradation  caused  by 
the  accumulation  of  cracks  in  a  straight  fibrous  composite,  we  proceed  in  the  next  section 
and  extend  this  procedure  to  the  case  of  undulated  fibers. 
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6.4  Damage  Analysis  in  Undulated  Fibrous  Composites 

The  stress  distribution  in  unidirectional  undulated  fiber  reinforced  composites  was 
found  in  chapters  4  and  5  for  a  variety  of  situations.  In  the  previous  section  the  procedure 
for  predicting  the  variation  of  the  composite  stiffness  with  crack  accumulation  in  straight 
fibrous  composites  was  introduced.  In  this  section  this  procedure  is  further  extended  to 
investigate  the  effects  of  damage  evolution  due  to  transverse  fiber  break  or  matrix 
cracking  in  undulated  fibrous  composites.  Following  the  same  procedure  of  chapter  4,  if 
we  start  by  an  undamaged  composite  segment  having  a  certain  fiber  geometric 
undulation,  we  identify  and  anal5n;ically  describe  local  tangents  to  the  fiber.  Second,  for  a 
given  applied  load  we  transform  the  global  coordinates  and  applied  loads  to  local 
coordinate  systems  consisting  of  the  tangent  directions  and  their  normals  in  the  plane  of 
undulation.  Next,  we  adapt  the  results  obtained  for  straight  fibers  to  "straight"  segmented 
fibers  along  the  tangents.  The  local  stresses  obtained  under  this  procedure  are  then 
supplemented  with  the  stress  contributions  that  inherently  rise  in  any  oriented  direction 
with  respect  to  the  loading  direction.  Once  the  state  of  stress  at  each  point  of  the  two 
media  is  found,  the  maximum  stress  failure  criterion  is  used  to  check  for  failure.  Once  a 
failure  takes  place  at  some  point,  the  stresses  in  the  new  composite  unit  cell  are  obtained 
again  under  the  new  boundary  conditions  (either  a  free  fiber  surface  or  a  free  matrix 
surface).  These  are  then  used  to  predict  if  new  fiber  breaks/  matrix  cracks  will  develop 
and  the  location  of  these  damage  surfaces.  The  process  is  repeated  for  an  increasing 
applied  load  and  the  break  density  values  are  predicted.  At  every  damaged  state,  the  total 
axial  strain  in  the  composite  is  evaluated  from  the  strain  in  each  segment  and  the 
reduction  in  the  longitudinal  stiffness,  as  expressed  by  Young’s  modulus,  is  determined. 
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6.4.1  Damage  Evolution 

As  the  composite  is  subjected  to  increasing  loads  at  its  ends,  damage,  in  terms  of 
number  of  fiber  breaks,  is  assumed  to  propagate  at  discrete  steps.  At  each  step,  the  stress 
distribution  between  adjacent  breaks  is  monitored  and  the  point  of  possible  future  break 
is  predicted  using  the  maximum  stress  criterion,  namely 

Zy  -  Max[a^j^ ]  =  0 ,  (6.5) 

where  ^  is  the  axial  strength  of  the  fiber  material. 

Assume  we  start  with  a  damage  free  composite  with  its  length  chosen  such  that  it 
spans  a  complete  number  of  fiber  waves,  say  four  for  example.  Because  both  the 
geometry  and  loading  are  symmetric,  the  point  of  maximum  stress  will  be  exactly  at  the 
middle  of  the  composite.  Upon  increasing  the  applied  load,  a  single  break  occurs  at  this 
point,  redistributing  the  stresses  everywhere.  No  further  breaks  take  place  until  the  load 
increases  enough  for  the  failure  criterion  to  be  satisfied  at  the  one  quarter  points  of  the 
composite,  resulting  again  in  symmetric  break  pattern.  The  last  mid-span  break  takes 
place  in  each  single  wave  as  the  load  is  increased  further.  Once  the  composite  is  divided 
into  half  waves,  which  are  not  geometrically  symmetric  with  respect  to  the  loading,  any 
further  breakage  will  not  necessarily  occur  at  the  middle  point  between  already  existing 
breaks.  Rather  it  will  depend  primarily  on  the  undulation  angle,  or  alternatively  on  the 
fiber  waviness  amplitude. 


6.4.2  Slustration  and  Discussion 


To  understand  how  damage,  in  terms  of  transverse  cracks  or  fiber  breaks 
accumulates,  we  consider  a  generic  unit  cell  made  of  the  same  fiber  and  matrix  materials 
of  the  examples  in  chapters  5.  The  unit  cell  is  initially  assumed  damage  free  and 
subjected  to  a  very  small  axial  loading.  Upon  increasing  of  the  load,  a  single  crack  takes 
place  as  shown  in  figure  6.4.  Due  to  redistribution  of  loading,  no  additional  cracks  occur, 
until  the  loading  value  increases  to  14.7  MPa,  where  4  cracks  are  predicted  at  very  small 
load  intervals.  This  is  in  contrast  to  the  straight  fibrous  case  where  only  one  crack  would 
occur  in  a  half  unit  cell.  We  also  notice  that  due  to  the  shape  of  the  stress  distribution, 
new  cracks  do  not  occur  close  to  the  middle  points  between  already  existing  cracks. 

As  another  illustration  of  the  procedure  capabilities  we  consider  an  undulated 
fibrous  composite  with  the  unit  cell  having  a  length  to  amplitude  ratio  of  10.  The 
mechanical  properties  of  the  E-glass  fiber  and  Epoxy  matrix  are  the  same  of  the 
illustration  of  section  6.3. 

Figure  6.5  shows  the  load-damage  curve  for  this  composite  system.  Initially  a 
single  transverse  crack  takes  place  in  the  weak  epoxy  material.  Upon  increase  of  loading 
and  monitoring  stresses,  other  cracks  at  several  locations  are  predicted.  It  is  important  to 
note  that  due  to  undulation  new  cracks  do  not  necessarily  occur  at  the  middle  point 
between  two  already  existing  cracks  as  was  the  case  in  the  straight  fibrous  case.  Results 
for  the  straight  fibrous  case  are  also  plotted  on  the  same  figure  for  comparison.  The 
degradation  of  the  composite  as  expressed  in  terms  of  reduction  in  Young’s  modulus  is 
demonstrated  in  figure  6.6  for  both  the  straight  fibrous  case  as  well  as  for  the  undulated 
case.  Both  values  of  Young’s  moduli  are  normalized  in  each  case  with  respect  to  its  value 
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Figure  6.5:  Evolution  of  transverse  cracks  in  a  generic  unit 
cell  of  undulated  fibrous  composite. 
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in  the  absence  of  damage.  The  figures  are  not  smooth  as  they  represent  the  damage 
accumulation  from  one  discrete  loading  step  to  another.  Although  the  value  of  Eo  in 
undulated  fibrous  composites  is  already  smaller  than  that  of  the  straight  case,  the  rate  of 
degradation  is  not  as  high.  For  example  at  a  crack  density  of  1  mm'\  the  undulated 
fibrous  composite  retains  30%  of  its  longitudinal  stiffness  as  compared  to  20%  in  the 
straight  fibrous  case.  Therefore,  in  applications  where  both  types  can  be  used  undulated 
composites  will  not  deteriorate  as  drastically  upon  damage  accumulation. 


Figure  6.6:  Comparison  between  the  load-damage  curves  for  an  E-glass/Epoxy  composite 
for  both  the  straight  fibrous  case  and  the  undulated  case. 
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Figure  6.7:  Comparison  between  the  degradation  curves  for  an  E-glass/Epoxy 
composite  for  both  the  straight  fibrous  case  and  the  undulated  case. 
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CHAPTER? 


APPLICATION:  VIBRATION  CHARACTERISTICS  OF  DAMAGED 
FIBROUS  COMPOSITES 

7.1  Introduction 

In  this  chapter,  and  as  a  further  illustration  of  the  utility  of  our  modeling, 
procedure,  we  turn  our  attention  to  the  study  of  vibration  characteristics  of  a  fibrous 
composite  panel.  This  modeling,  which  was  developed  in  chapter  2,  is  general  and  can  be 
used  in  solving  many  classes  of  problems,  including  wave  propagation  and  vibrations. 
The  former  problem  received  a  continuous  attention  in  research,  see  for  example  Nayfeh 
[1],  Nayfeh  and  Nagy  [31],  Simmons  et  al.  [32]  and  Thurston  [33].  In  those  papers,  both 
approximate  and  exact  formulations  were  presented  for  the  interaction  of  fibrous 
composites  with  a  propagating  wave.  Much  less  work  is  available  on  the  influence  of 
microstructure  on  the  vibration  characteristics  of  fibrous  composites,  even  in  the  absence 
of  damage.  Here,  we  extend  our  micromechanical  model  in  order  to  investigate  these 
characteristics  with  and  without  damage.  The  basic  three  possible  fiber/matrix 
arrangements  are  shown  in  figure  7. 1 .  Damage  will  be  simulated  by  either  a  fiber  break,  a 
matrix  crack  or  by  combinations  at  selective  locations  along  the  fiber  direction.  Some 
anal)^ical  expressions  for  the  resonance  frequency  modes  will  be  derived.  Numerical 
results  will  be  presented  mostly  as  comparisons  between  damaged  and  undamaged 
situations. 


^4 


12  Mathematical  Formulation 


Since  this  vibration  application  is  intended  only  for  illustration,  and  to  avoid 
unnecessary  complicated  equations,  we  will  make  an  extra  simplifying  assumption  to  the 
micromechanical  model.  We  will  neglect  the  small  effect  of  the  dvidz  on  the  shear 
constitutive  relation,  Eq.  (2.9)  for  both  the  fiber  and  matrix  materials,  compared  to  the 
significant  direct  effect  of  du  I  dr .  This  decouples  the  integrated  forms  of  this  equation, 
namely  equations  (2.30a)  and  (2.30,b),  into 

and 

where 


G  = 


8/55^55 

''2^(»*55  -/550 


(7.3) 


Now,  specializing  the  constitutive  relations  (2.38a)  and  (2.38b)  to  r  =  r,  and  invoking  the 
interface  continuity  equation  (2.10)  yields 


IrifD  dz 


(7.4) 


where 


m 
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Combining  this  equation  with  (7.2),  we  solve  for  B  in  terms  of  the  strain  components 


dUf 

and 

dz 


dz 


as 


5  = 


In. 


m. 


'55 


du^ 

dz 


-fssQ- 


dUf 

dz 


where 


(7.6) 


^("*55  -/ssG)  ■ 


(7.7) 


Finally,  substituting  from  this  equation  back  into  the  averaged  constitutive 
relations  (2.26)  and  (2.27),  combining  the  resulting  equations  with  equations  (2.24)  and 
(2.25)  and  using  equation  (7.1)  to  eliminate  A,  leads  to  the  two  coupled  partial  differential 
equations; 


oz 


(7.8) 


-G(m„-m^)  =  0. 


(7.9) 


Here  the  coefficients  Cj,,  Cjj,  Cjj,  Cjj  are  given  as  combinations  of  the  fiber  and  matrix 
properties  and  volume  fractions  as 

<^11  =  W//11  -  fufssFQ ,  ^12  =  .  (7.10a) 

C21  =  ~^xifiiFQ ,  C22  =  ~  ^u^ssF  ■  (7.10b) 

The  behavior  of  the  composite  is  completely  described  by  equations  (7.8)  and 

(7.9).  As  seen,  this  system  consists  of  two  partial  differential  equations  that  couple  the 
two  longitudinal  averaged  displacements.  These  equations  retain  the  integrity  of  the 
distribution  of  displacements  and  stresses  in  the  fiber  and  matrix  individually  subject  to 
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the  influence  of  the  interfacial  transfer  couplings.  Thus,  upon  loading  the  composite,  the 
distribution  of  stresses  and  displacements  in  the  fiber  and  matrix  components  are  readily 
calculable.  Furthermore,  the  form  of  the  coupled  equations  in  our  model  makes  it 
possible  to  incorporate  damage  mechanisms  such  as  fiber  break  and  matrix  crack  in  the 
composite  system.  This  system  is  general  and  can  be  applied  to  a  variety  of  physical 
situations. 

7.3  Applications 

In  the  remainder  of  this  chapter  we  apply  the  present  modeling  to  the  study  of  the 
vibration  characteristics  of  fibrous  composite  systems  with  or  without  damage.  For 
damage  fi'ee  composite  systems,  we  choose  a  generic  composite  panel  that  occupies  the 
region  -Lll<z<Lll  with  its  center  located  at  z=0.  Damage  is  then  simulated  as  a 
stress  free  boundary  in  either  the  fiber  (for  fiber  break)  or  the  matrix  (for  matrix  crack) 
along  the  panel  length. 

Depending  upon  the  number  of  damage  zones,  the  panel  is  divided  into  discrete 
sub-panels  connected  at  the  damage  locations.  Thus,  for  an  arbitrary  number  of  damage 
locations,  A'^+l  sub-panel  segments  are  required.  Each  segment  has  the  individual  length 
Z-/,  /=1,  2,  ...iV+1,  totaling  L.  Furthermore,  we  assign  for  each  segment  /  the  local 
coordinate  Zi  with  origin  located  at  its  center;  thus  it  extends  from  -L-jl  to  Z,;/2.  Since 
equations  (7.8)  and  (7.9)  hold  for  all  segments,  formal  solutions  for  each  can  then  be 
written  in  terms  of  two  unknown  amplitudes.  By  satisfying  the  appropriate  continuity 
conditions  across  segment  interfaces  and  followed  by  invoking  the  stress  fi’ee  conditions 
at  the  left  and  right  ends  of  the  total  panel,  we  obtain  the  characteristic  equation  for  the 


‘78 


resonance  frequency.  Resonance  frequencies  for  situations  involving  other  than  stress- 
free  ends  of  the  panel  can  be  just  as  easily  obtained  and  will  not  be  presented  in  this 
paper. 

Similar  to  what  we  did  in  all  cases  involving  dynamic  loading,  a  formal  solution 
of  equations  (7.8)  and  (7.9)  can  be  assumed  in  the  time  harmonic  form 

(7.11a) 

(7.11b) 

This  results  in  a  coupled  system  of  equations  in  the  amplitudes  Xi  and  X2,  rearranged  in 
matrix  form  as 

n^Pj-ct)^  +  Cj,a^  -  G  +  G 

-  +  G  n„p„a)^  +  -  G 

For  nontrivial  solutions,  the  determinant  of  equation  (7.12)  must  be  zero.  For  a  given 

frequency,  the  resulting  characteristic  equation  admits  four  solutions  for  a,  or  more 

precisely  two  solutions  for  a^.  Thus,  solutions  for  a  occur  in  two  pairs,  each  pair 

having  two  a ’s  that  are  negative  of  each  other.  We  shall  label  these  roots  as  a^, 

2,  3,  4,  and  further  adopt  the  convention  that  cZj  “  ~  > 

equation  (7. 12)  can  be  used  to  obtain  the  amplitude  ratio 

(713) 

Using  the  method  of  superposition,  the  formal  solutions  for  the  displacements  and 
stresses  become 


lx. 


=  0. 


(7.12) 
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(7,15) 


From  the  definitions  of  V2q,  Dig,  D2q,  it  is  clear  that 


^3+1  ~^q’  ^\q*\  ~  ^\q  ^2q+l  ~  ^2q’  ^1>3. 


(7.16a) 

(7.16b) 

(7.16c) 

(7.17) 


The  formal  solutions  (7.14)  and  (7.15)  are  complete  when  the  four  values  of  Xig 
are  determined.  Four  displacement  and/or  stress  boundary  conditions,  depending  upon  the 
problem  at  hand,  are  therefore  required.  For  free  longitudinal  vibration  of  an  undamaged 


panel,  for  example,  we  require  the  vanishing  of  the  normal  stresses  at  both  ends  z  =  L  /  2 


and  z  =  -L  !  2. 


7.3.1  Undamaged  Fiber-Matrix  Composite 

In  the  absence  of  damage,  figure  7.  la  shows  that  only  one  segment  is  necessary  to 
model  the  composite  panel.  Imposing  the  stress  free  conditions,  namely  5^  =  0  and 

at  z  =  -LI2  and  z  =  Z,/2,  on  the  formal  solution  (7.15),  leads  to  the 
characteristic  equation 

(e“.i  -  )(e“^  -  )  =  0 .  (7.18) 
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7.3.2  Case  of  a  Single  Fiber  Break 

For  cases  involving  a  single  fiber  break,  say  at  z  =  ,  as  shown  in  figure  7.1b, 

then,  according  to  our  convention,  the  panel  is  composed  of  two  segments.  The  formal 
solutions  (7.14)  and  (7.15)  are  then  adapted  for  each  segment.  At  z,  =4/2  and 
Zj  =  -Zj  12,  continuity  of  the  longitudinal  displacement  and  stress  are  required, 
whereas  the  stress  vanishes  for  both  segments  at  that  location.  At  the  same  time, 
imposing  the  stress  free  conditions  at  the  outer  ends  of  the  system,  namely  -  0  and 


=0  at  Zj  =  -Z,  /2  and  at  results,  after  a  set  of  row  operations  on  the 

resulting  determinant  of  the  boundary  conditions,  in  the  characteristic  equation 


f^sAiC 


C(l)  C(2)  0(1)  0(2) 


(7.19) 


where 


gU)  _ 


/=1,3  ,  j=\,2. 


7.3.3  Case  of  a  Single  Matrix  Crack 

If  instead  of  a  fiber  break  we  have  a  matrix  crack  at  the  location  z  =  Zj ,  as 
demonstrated  in  figure  7.1c,  then  the  interface  conditions  at  Zj  =  Z,  /2  and  =  -Zj  /2 
require  the  continuity  of  the  fiber  displacement  Wy  and  the  fiber  stress  and  the 

vanishing  of  the  matrix  stress  in  both  segments  at  the  crack  location.  By  inspection, 
we  derive  the  characteristic  equation  for  the  case  of  matrix  crack  by  simply  replacing  Du 
by  D2i  and  eliminating  V-,,  /=1,3,  in  equation  (7. 19)  as 
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(7.20) 


C(l)  C(2) 

-aX.  ‘J'S 


D  (  ^  ^-D  {  ^  e'°‘^  + 

^2Acr(2)‘^  ^  r-fn ^  -^zaVom*^  ^ 


C(l)  C(2) 


)  =  0. 


7.3.4  Case  of  a  Damage  Combination 

The  above  procedure  can  be  used  to  obtain  solution  for  the  cases  involving 
multiple  damage  in  a  composite  panel.  In  a  case  involving  n  segments,  connected  at 
damage  locations,  a  4n  x4w  determinant  need  to  be  manipulated.  For  n  larger  than  2,  the 
algebraic  manipulation  of  the  determinant  becomes  increasingly  tedious.  Hence,  we 
resort  to  numerical  methods  to  find  the  zeros  of  the  An  x  An  determinant  reflecting  the 
appropriate  interface  and  boundary  conditions. 

7.4  Discussion  and  Numerical  Illustrations 

Casual  inspection  of  the  solution  of  the  undamaged  composite,  equation  (7.18), 
reveals  that  two  decoupled  sets  of  vibration  harmonics  are  possible  to  exist  in  the  present 
panel.  This,  however,  requires  both  values  of  a  to  be  pure  imaginary.  In  real  situations, 
this  is  not  always  the  case.  We  find  that,  for  sufficiently  low  resonance  frequencies,  only 
one  value  of  a  is  imaginary,  reflecting  a  vibrating  mode,  while  the  other  value  is  real, 
reflecting  an  evanescent  mode.  Whether  a  is  real  or  imaginary  depends  on  the  resonance 
frequency  ranges.  Our  discussion  of  this  point,  and  of  others  to  follow,  can  best  be 
facilitated  with  the  aid  of  a  representative  numerical  illustration. 

In  all  of  our  calculations,  we  shall  choose  graphite-epoxy  as  a  representative 
composite  panel.  This  panel  has  the  following  necessary  geometric  and  material 
parameters.  The  total  length  and  radius  of  the  fiber  are  arbitrarily  chosen  as  20  cm  and  1 
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cm,  respectively  with  a  fiber  volume  fraction  of  0.55.  The  properties  of  both  the  graphite 
and  the  epoxy  components  are  quoted  from  reference  [1]  and  listed  for  convenience  in 
stiffness  matrix  format  as: 


235.0  3.69  3.69 
26.0  3.32 
26.0 


0 

0 

0 

5.52 


0 

0 

0 

0 

28.2 


0 

0 

0 

0 

0 

28.2 


'sym 


xl0‘° 


2 

dyne/cm 


for  graphite,  and 


m..  = 

V 


7.67  3.75  3.75  0  0 

7.67  3.75  0  0 

7.67  0  0 

1.96  0 

1.96 


0  ■ 
0 
0 
0 
0 

1.96 


\syni 


X  10'°dyne/cm^ 


for  epoxy.  Also  needed  are  the  material  densities  Pf  =  1 .79  gm/cm^  for  graphite  and 

=  1.26  gm/cm  for  epoxy.  To  further  facilitate  our  discussion,  we  remind  the  reader 
of  the  trivial  solutions  for  or  in  a  single  medium.  For  single  material  panels,  a  is 
calculated  simply  from  a  -  /  p ,  where  Cn  and  p  are  the  longitudinal  stiffness  and 

density,  respectively. 

While  for  the  simple  case  of  the  single  medium  a  is  linear  in  co,  this  is  not  true  in 
the  composite  problem  and,  in  fact,  equation  (7.12)  reveals  that  a  has  a  quadratic 
dependence  on  u);  this  is  illustrated  numerically  in  figure  7.2.  Also  included  for 
comparison  in  this  figure  are  the  a ’s  of  the  single  graphite,  single  epoxy  and  mixture 
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media  and  designated  as  and  _  respectively.  The  mixture  value  is  based 

upon  the  effective  longitudinal  property  (along  the  fiber  direction)  of  the  composite.  This 
property  is  obtained  from  the  original  mixture  system  in  the  static  limit.  In  this  limit,  the 
average  displacements  in  both  the  fiber  and  matrix  approach  each  other.  By  summing 
both  equations  subject  to  i/y  =  =  m  ,  we  obtain  the  single  wave  equation 


d^u 

(Ci,  +Ci2  +C2,  +C22)^  =  (nfPf  +nn,Pn,) 


dz^  V  /  -  y  m,  ’ 


(7.21) 


reflecting  the  homogenized  medium  with  the  effective  stiffness  and  density  given 
respectively  as 


Figure  7.2:  Variation  of  a  with  the  radial  frequency  a;  for  undamaged  panel. 


^4 


(7.22a) 


Thus  ccmix  is  calculated  from 

a  =  • 


(7.22b) 


(7.23) 


Inspection  of  figure  7.2  reveals  that,  at  relatively  low  frequency  ranges,  only  one 
of  the  composite  a ’s  is  pure  imaginary  whereas  the  other  is  pure  real  signifying  the 
existence  of  a  single  set  of  vibrating  harmonics  and  that  their  values  are  relatively  close 
to  that  of  the  effective  mixture  values.  Here,  the  composite  behaves  as  a  single  effective 
(homogenized)  material  with  both  of  its  constituents  loose  their  individual  identity  and 
move  in  concert  with  each  other.  As  the  frequency  increases  beyond  a  critical  value,  both 
a ’s  become  pure  imaginary  signifying  the  existence  of  two  sets  of  vibrating  modes.  The 
analytical  expression  for  this  critical  frequency  is  obtained  by  setting  a  =  0  in  the 
determinant  of  equation  (7. 12),  and  solving  for  o,  leading  to: 


p(»/P/+”m~Aj 


(7.24) 


For  the  present  system,  this  critical  frequency  is  calculated  as  6.7274  x  10^  rad/s. 

The  evolution  of  the  two  sets  of  harmonics  for  increasing  frequency  is  interesting. 
For  moderate  frequency  ranges  beyond  the  critical  frequency,  these  two  modes  stay 
interdependent  in  that  each  is  dependent  on  both  fiber  and  matrix  properties  and  the 
composite  geometry.  As  the  frequency  increases  further,  the  two  sets  become  more  and 
more  distinct  and,  as  the  frequency  approaches  infinity,  they  become  identifiable  with 
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those  of  the  pure  fiber  and  pure  matrix  panel  modes.  This  can  also  be  seen  from  figure 
7.2  where  the  two  branches  asymptotically  approach  those  of  and  .  In  the 

meantime  is  seen  to  diverge  a  little  more  from  the  lower  branch  suggesting  that  the 
mixture  mode  is  only  valid  in  the  low  frequency  ranges. 

It  can  be  seen  that  the  high  frequency  limits  also  have  their  roots  in  the  original 
mixture  equations  (7.8)  and  (7.9).  As  the  frequency  is  increased  beyond  limit,  the 
characteristic  equation  (7.18)  uncouples  and  leads  to  two  limiting  values  for  a ’s.  Except 
for  a  slight  coupling  from  the  and  C21  terms,  these  two  limiting  values  coincide  with 
Uf  and  a„ . 

Having  given  a  qualitative  discussion  of  the  behavior  of  the  vibrating  modes,  we 
now  move  to  present  numerical  calculations  for  cases  in  which  damage  is  either  absent  or 
present.  We  start  with  a  damage  free  panel.  For  stress  free  outer  ends,  the  resonance 
frequencies  are  obtained  from  satisfying  equation  (7. 1 8).  For  this  specific  situation,  these 
frequencies  can  also  be  extracted  from  figure  7.2.  It  should  be  noted  that  this  figure  is 
scaled  such  that  the  resonance  frequencies  are  obtained  directly  by  projecting  the 
intersections  of  the  horizontal  lines  and  the  a  curves  on  the  abscissa.  The  lowest 
(fundamental)  resonance  frequency  co^  for  this  undamaged  case  can  either  be  roughly 

estimated  from  the  figure  or  exactly  calculated  from  equation  (7. 18)  as  1 .4355  x  10^  rad/s. 
Higher  resonance  frequencies  of  this  panel  can  also  be  easily  obtained.  Here  we  list  the 
first  fifteen  and  label  them  consecutively  as  in  Table  7.1.  For  future 

reference,  we  list  their  normalized  values  with  respect  to  the  fundamental  frequency  cUj. 
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Break  Location  (cm) 

Figure  7.3;  Variation  of  the  normalized  natural  frequencies  coid,  Qhd,  coid  and  cond 
with  the  damage  location  for  the  case  of  broken  fiber. 

resonance  frequencies.  This  is  due  to  the  fact  that  the  graphite  fiber  is  much  stiffer  than 
the  epoxy  and  thus  little  influence  can  be  attributed  to  the  cracked  matrix. 

As  an  illustration  of  the  versatility  of  the  present  analysis,  we  finally  investigate  the 
sensitivity  of  the  resonance  frequencies  to  the  relative  volume  fractions  of  the  fiber  and 
matrix  constituents.  The  variation  of  the  first  three  fundamental  fi'equencies  of  the 
composite  panel  with  the  fiber  volume  fraction  is  depicted  in  figure  7.4  in  the  presence 
and  absence  of  damage.  For  this  specific  illustration,  damage  is  taken  in  the  form  of  a 
fiber  break  at  the  center  of  the  20  cm  long  and  1  cm  radius  fiber  composite  panel.  The 
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dotted  and  solid  curves  belong  to  damage  free  and  damaged  situations.  Inspection  of 
figure  7.4  leads  to  some  interesting  conclusions.  First  the  two  limits  of  rij^  =Q  and  1 
belong  to  a  panel  made  up  of  only  matrix  or  fiber  material,  respectively.  At  the  limit  of 
zero  fiber  volume  fraction,  the  fundamental  frequencies  of  the  damaged  and  undamaged 
panels  coincide. 


Figure  7.4;  Fiber  volume  fraction  dependence  of  the  first  three  natural  frequencies  of  a 
composite  panel  for  the  undamaged  case  and  for  the  case  of  a  fiber  break  at  z=0. 


As  the  fiber  volume  fraction  increases,  emphasizing  thinner  matrix  cladding,  the 
solutions  for  the  damaged  panel  depart  away  from  the  corresponding  undamaged  one.  In 


the  limit  as  the  fiber  volume  fraction  approaches  unity,  the  fundamental  mode  approaches 
zero,  whereas  the  third  mode  approaches  that  of  the  second.  This  interesting  limit  can  be 
easily  explained  by  investigating  the  behavior  of  the  corresponding  mode  shapes.  For  the 
sake  of  illustration  we  investigate  the  evolution  of  the  third  mode  shape  with  varying 

fiber  volume  fraction.  As  can  be  seen  from  figure  7.5,  when  the  fiber  volume  fraction 
increases,  the  two  nodal  points  close  to  the  center  approach  the  crack  location  at  mid- 
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Figure  7.5:  Variation  of  the  third  matrix  stress  mode  shape  with  the  fiber  volume 
fraction  for  the  case  of  a  broken  fiber  atz=10. 
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span.  Meanwhile,  the  matrix  stress  peak  at  the  crack  location  grows,  a  behavior  consistent 
with  the  fact  that  the  ever-thinner  matrix  still  has  to  carry  the  same  load  across  the  crack. 
At  the  extreme  limit  =  1 ,  this  evolution  results  in  the  third  mode  being  coincident  with 
the  second  mode. 

I 

In  this  chapter,  the  newly  developed  micromechanical  model  is  applied  to  the  study  the 
vibration  characteristics  of  fibrous  composites  with  damage.  The  model  proved  to  be 
easily  adapted  for  cases  involving  single  fiber  break,  single  matrix  crack  or  multi  damage 
situations.  It  is  used  to  investigate  the  dependence  of  the  composite  panel  resonance 
frequencies  and  mode  shapes  on  the  damage  location  and  on  the  fiber  volume  fraction  in 
the  case  of  a  single  fiber  break.  Results,  which  are  also  presented  in  [34],  show  several 
interesting  characteristics  as  the  fiber  volume  fi-action  approaches  its  limiting  values. 
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CHAPTER  8 


CONCLUSION 

In  this  work,  a  unified  micromechanical  model  for  the  prediction  of  the  stress 
distribution  in  undulated  fibrous  composites  was  developed  and  used  to  study  damage 
accumulation  in  this  important  class  of  composites.  First,  the  unidirectionally  undulated 
fibrous  unit  cell  is  divided  into  segments  each  with  a  local  principal  direction  parallel  to 
the  fiber  tangent  in  this  segment.  Then  the  interfacial  values  of  some  of  the  stress  and 
displacement  components  were  used  together  with  the  averaging  technique  to  reduce  the 
system  of  resulting  equations  into  a  quasi-one-dimensional  one.  In  doing  so,  certain 
approximations  had  to  be  made  regarding  the  dependence  of  some  of  the  stress  and 
displacement  components  on  the  radial  direction.  The  resulting  fourth  order  differential 
equation  was  solved  for  the  stresses  in  the  “local”  fiber  direction,  and  the  resulting  shear 
stress  was  then  enhanced  with  the  already  existing  stress  between  the  fiber  and  the 
matrix.  Damage  in  the  composite  in  the  form  of  a  fiber  break  or  a  matrix  crack  is 
modeled  by  the  equation  through  using  the  right  set  of  boundary  conditions.  The  accuracy 
of  the  modeling  is  verified  by  several  comparisons  with  both  available  published  results 
and  finite  element  generated  solutions. 

Having  obtained  the  stress  distribution  in  both  constituents,  the  maximum  stress 
failure  criterion  is  used  to  check  if  additional  damage  will  occur  in  any  of  the  constituents 
upon  the  increase  of  loading.  The  process  is  repeated  for  as  many  load  steps  until 
catastrophic  failure  occurs.  At  each  load-step,  the  integration  of  the  displacements  in  the 
global  composite  axial  direction  provides  a  means  of  estimating  the  composite  stiffness. 
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Monitoring  stiffness  and  damage  density  upon  loading  the  composite  provides  important 
performance  information.  This  includes  the  load-damage  diagram  and  the  stiffness- 
damage  diagram.  These  diagrams  stress  the  fact  that  weakening  due  to  damage 
accumulation  is  more  significant  in  composites  reinforced  with  straight  fibers  than  in 
those  with  undulated  ones 

A  final  application  of  the  model  was  that  in  the  area  of  vibrations.  The  model  was 
easily  adapted  to  study  vibrations  of  damaged  straight  fibrous  composites.  The  fact  that 
damage  location  affects  the  system  natural  frequencies  and  mode  shapes  is  of  practical 
significance,  since  it  can  be  used  in  the  prediction  of  damage  location. 

The  proposed  formulation  opens  several  paths  for  future  studies.  From  geometry 
point  of  view,  the  analysis  can  be  extended  to  cover  both  the  warp  and  weft  yams. 
Existing  models  which  deals  with  both  yams  tend  to  neglect  the  continuity  of  their 
undulation,  as  for  example,  the  mosaic  model.  Although  not  discussed  here,  the  effect  of 
undulation  on  the  degradation  in  transverse  properties  can  also  be  investigated.  Finally, 
other  potential  capabilities  for  the  model  need  to  be  fiilly  explored.  This  includes 
problems  of  wave  propagation  and  vibrations  in  textiles. 
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